February 1,2008 



o 
o 

(N 

< 

(N 



Oh 

6 



> 

00 
(N 

q 

00 

o 
i> 
o 



1— 1 

>< 



I 



Pulsation tomography of rapidly oscillating Ap stars* 
Resolving the third dimension in peculiar pulsating stellar atmospheres 

T. Ryabchikova 1 ' , M. Sachkov 2 , 0. Kochukhov 3 , and D. Lyashko 4 

Department of Astronomy, University of Vienna, Tiirkenschanzstrasse 17, A-1180 Wien, Austria 
Institute of Astronomy, Russian Academy of Sciences, Pyatnitskaya 48, 1 19017 Moscow, Russia 
Department of Astronomy and Space Physics, Uppsala University Box 515, SE-751 20 Uppsala, Sweden 
Tavrian National University, Yaltinskaya 4, 330000 Simferopol, Ukraine 



Received / Accepted 

ABSTRACT 

Aims. We present detailed analysis of the vertical pulsation mode cross-section in ten rapidly oscillating Ap (roAp) stars based on 
spectroscopic time-series observations. The aim of this analysis is to derive from observations a complete picture of how the amplitude 
and phase of magnetoacoustic waves depend on depth. 

Methods. We use the unique properties of roAp stars, in particular chemical stratification, to resolve the vertical structure of p-modes. 
Our approach consists of characterising pulsational behaviour of a carefully chosen, but extensive sample of spectral lines. We analyse 
the resulting amplitude-phase diagrams and interpret observations in terms of pulsation wave propagation. 

Results. We find common features in the pulsational behaviour of roAp stars. Within a sample of representative elements the lowest 
amplitudes are detected for Eun (and Fe in 33 Lib and in HD 19918), then pulsations go through the layers where Ho- core, Nd, 
and Pr lines are formed. There RV amplitude reaches its maximum, and after that decreases in most stars. The maximum RV of the 
second REE ions is always delayed relative to the first ions. The largest phase shifts are detected in Tbm and Thin lines. Pulsational 
variability of the Thin lines is detected here for the first time. The Yn lines deviate from this picture, showing even lower amplitudes 
than Eun lines but half a period phase shift relative to other weakly pulsating lines. We measured an extra broadening, equivalent to 
a macroturbulent velocity from 4 to 1 1-12 km s -1 (where maximum values are observed for Tbm and Thm lines), for pulsating REE 
lines. The surface magnetic field strength is derived for the first time for three roAp stars: HD 9289 (2 kG), HD 12932 (1.7 kG), and 
HD 19918 (1.6 kG). 

Conclusions. The roAp stars exhibit similarity in the depth-dependence of pulsation phase and amplitude, indicating similar chemical 
stratification and comparable vertical mode cross-sections. In general, pulsations waves are represented by a superposition of the 
running and standing wave components. In the atmospheres of roAp stars with the pulsation frequency below the acoustic cut-off 
frequency, pulsations have a standing-wave character in the deeper layers and behave like a running wave in the outer layers. Cooler 
roAp stars develop a running wave higher in the atmosphere. In stars with pulsation frequency close to the acoustic cut-off one, 
pulsation waves have a running character starting from deep layers. The transition from standing to running wave is accompanied by 
an increase in the turbulent broadening of spectral lines. 

Key words, stars: atmospheres - stars: chemically peculiar - stars: magnetic fields - stars: oscillations 

1. Introduction More than 30 cool Ap stars exhibit high-overtone, low- 

degree, non-radial p-mode pulsations with periods in the range 

About 10 % to 20 % of upper main sequence stars are charac- of 6 _ 21 minutef , (Kurtz & Martinez 2QQOu with their observe d 
tensed by remarkably rich line spectra, often containing numer- lsation amplitudes modulated according to the visible mag- 

£ . ous unidentified features. Compared to the solar case, overabun- netic fieM stmcture xhese m osci u ating Ap (roAp) stars 

dances of up to a few dex are often inferred for some iron-peak m k objects for asteroseismologVi which pre sently is the most 

and rare-earth elements, whereas some other chemical elements powerful tool for testing theories of stellar stmcture and evolu . 

tion. 



are found to be underabundant (Ryabchikova et al. 2004). Some 



of these chemically peculiar (Ap) stars also exhibit organised Recent jn observational studies of roA stars has 

magnetic fields with a typical strength of a few kG. Chemical been achieyed fe consideri hi h time - r esolution spectroscopy 

peculiarities are believed to result from the influence of the mag- fa a(Mition tQ ^ hi h d photometric mea sure- 

netic field on the diffusing ions, possibly in combination with mentf , Hi h alit time . resolve d measurements of magnetic 

the influence of a weak, magnetically-directed wind (e.g., Babel pulsatorfj haye reyealed a surprising diversity in the pulsational 

' "■ behaviour of different lines in the roAp spectra (e.g., Kanaan & 

Send offprint requests to: T. Ryabchikova, HatzesHMl Kurtz et al. [2006]and references therein). Detailed 

e-mail: ryabchik@inasan.ru analyses of the bright roAp star y Equ (Savanov et al. 119991 

* Based on observations made with the SAO 6-m telescope, with the Kochukhov & Ryabchikova 2001a) demonstrate that spectro- 

Canada-France-Hawaii Telescope, and the ESO VLT (DDT programme scopic pulsational variability is dominated by the lines of rare- 

274.D-5011 and programme 072.D-0138, retrieved through the ESO earth ions, especially those of Pr and Nd, which are strong 

archive). and numerous in the roAp spectra. On the other hand, light 



and iron-peak elements do not pulsate with amplitudes above 
50-100 ms , which is at least an order of magnitude lower 
in comparison with the 1-5 kms~' variability observed in the 
lines of rare-earth elements (REE). Many other roAp stars were 
proven to show very similar overall pulsational behaviour (e.g., 
K ochukh ov & Ryab chikovaBOOlbl Balona 120021 Mkrtichian et 
al. 2003; Kurtz et al. 2005a), with an exceptional case in 33 Lib, 
which shows small-amplitude variations in Fe lines, and possibly 
two other stars, HD^ 12932 and HD 19918, where RV variations 
in the Fei 5434.52 A line have been marginally detected (Kurtz 
et al. l2005bb . 

The early spectroscopic studies speculated that the unusual 
diversity of the pulsation signatures in roAp spectra can be at- 
tributed to an interplay between the p-mode pulsation geometry 
and inhomogeneous horizontal or vertical element distributions 
(see discussions in Kanaan & Hatzes fT998l Savanov et al.[l999; 
Baldry & Bedding |2000l Mkrtichian et al. [2000l . However, none 
of these studies present models capable of explaining pulsations 
in real stars. Peculiar characteristics of the /9-mode pulsations in 
a roAp star were finally clarified by Ryabchikova et al. (f2002), 
who were the first to empirically determine vertical stratifica- 
tion of chemical elements and relate chemical profiles to pulsa- 
tional variability. In their study of the atmospheric properties of 
y Equ, Ryabchikova et al. show that the light and iron-peak el- 
ements are enhanced in the lower atmospheric layers, whereas 
REE ions are concentrated in a cloud with a lower boundary at 
log rgooo S= -4 (Mashonkina et al. 2005). Thus, high-amplitude 
pulsations observed in REE lines occur in the upper atmosphere, 
while lines of elements showing no significant variability form in 
the lower atmosphere. This leads to the following general picture 
of roAp pulsations: we observe a signature of a magnetoacous- 
tic wave, propagating outwards through the chemically stratified 
atmosphere with increasing amplitude. 

In addition to remarkable pulsational behaviour, the REE 
lines formed in the upper atmospheric layers of roAp stars ex- 
hibit an extra broadening, corresponding to a macroturbulent 
velocity Vm aC ro — 10 kms -1 , which cannot be attributed to the 
chemical stratification or magnetic field effects (Kochukhov & 
Ryabchikova 2001a). In the recent detailed line-profile variabil- 
ity study, Kochukhov et al. (|2007 ) have presented evidence for 
the existence of peculiar asymmetric oscillation patterns in the 
broad REE lines of several roAp stars. It was demonstrated that 
the inferred pulsation signatures can be reproduced with the 
spectrum synthesis calculations, which contain an extra pulsa- 
tional line width variability in addition to the usual velocity per- 
turbations. These results suggest that a turbulence zone modu- 
lated by pulsations probably exists in the upper atmospheres of 
roAp stars. 

The presence of significant phase shifts between pulsa- 
tion radial velocity (RV) curves of different REEs (Kochukhov 
& Ryabchikova 2001a) or even lines of the same element 
(Mkrtichian et al. 2003 ) can be attributed to the chemical stratifi- 
cation effects and, possibly, to a short vertical length of running 
magnetoacoustic wave. However, it is not immediately clear 
whether all spectroscopic observations of roAp stars can be fit- 
ted into this simple picture and to what extent magnetoacoustic 
pulsation theories can explain these observations. The wide di- 
versity of pulsation signatures (in particular, phases and bisector 
variability) of the REE lines probing similar atmospheric heights 
and, especially, the presence of pulsation node in the atmosphere 
of 33 Lib (Mkrtichian et al. 2003 ) are inexplicable in the frame- 
work of the non-adiabatic pulsation models developed by Saio & 
Gautschy (2004) and Saio (2005). These calculations correctly 
predict an increase in pulsation amplitude with height, but show 



neither nodes nor rapid phase changes in the REE line-forming 
region. To resolve this discrepancy it is therefore imperative to 
analyse in detail pulsational variations of many different ions in 
the spectra of representative sample of roAp stars. Only in this 
way can we derive meaningful observational constraints for pul- 
sation theories and search for regular patterns in the pulsation 
characteristics of different roAp stars. 

First general results for roAp stars were presented by Kurtz 
at al. d2005al) . They considered bisector behaviour of the Ha 
core and the Ndm6145 A line in 10 roAp stars. Based on the 
Ha core measurements, the authors point out that the increasing 
amplitude with height in the atmosphere is a common character- 
istic of all stars in their sample. Kurtz et al. (2005a) also note 
that in some stars, e.g. HD 12932, pulsations have a standing- 
wave behaviour, while more complex pulsations are observed in 
some other stars. These results were obtained from the analy- 
sis of bisector variability of one particular line, Ndm6145 A. 
This feature is not optimal for the bisector analysis because of 
the different blending effects both in the line core and the line 
wings that depend on the effective temperature and/or chemical 
anomalies. 

In this study we have embarked on the task of obtaining a 
detailed vertical cross-section of the roAp pulsation modes. We 
used the unique properties of roAp stars, in particular chemical 
stratification, to resolve the vertical structure of /?-modes and to 
study propagation of pulsation waves at a level of detail, that was 
previously only possible for the Sun. Our pulsation tomography 
approach consists of characterising the pulsational behaviour of 
a carefully chosen, but extensive, sample of spectral lines in- 
cluding weak ones and subsequently interpretating the obser- 
vations in terms of pulsation wave propagation. Therefore, our 
sample was limited to slowly rotating roAp stars. The aim of this 
analysis is to derive observationally a complete picture of how 
the amplitude and phase of magnetoacoustic waves depend on 
depth and to correlate resulting vertical mode cross-sections with 
other pulsational characteristics and with the fundamental stel- 
lar parameters. Furthermore, we envisage that unique 3-D maps 
of roAp atmospheres and non-radial pulsations can be derived 
by combining pulsation tomography results presented here with 
the horizontal pulsation maps obtained with Doppler imaging 
(Kochukhov 2004) or moment analysis (Kochukhov 2005). 

Our paper is organised as follows. Section [2] describes ac- 
quisition and reduction of roAp time-series spectra. The choice 
of targets and physical properties of the stellar sample are dis- 
cussed in Sect. [3] Radial velocity measurements and period anal- 
ysis are presented in Sect. [4] Pulsation tomography results are 
summarised in Sect. [6] Bisector variability is analysed in Sect. [5] 
Results of our study are summarised and discussed in Sect. [7] 

2. Observations and data reduction 

The main observational dataset analysed in our study consists of 
958 observations of 8 roAp stars, obtained with the UVES in- 
strument at the ESO VLT between October 8, 2003 and March 
12, 2004 in the context of the observing programme 072. D- 
0138 (Kurtz et al. [2006t . The ESO Archive facility was used 
to search and retrieve science exposures and the respective cal- 
ibration frames. Observations of each target covered 2 hours 
and consisted of an uninterrupted high-resolution spectroscopic 
time-series with a total number of exposures ranging from 69 
to 265. The length of individual exposures was 40 s or 80 s , ex- 
cept for the brightest roAp star HD 128898 (a Cir), for which a 
1.5 s exposure time was used. The ultra-fast (625kHz/4pt) read- 
out mode of the UVES CCDs allowed us to limit overhead to 



Table 1. Log of time-series observations for roAp stars. 



Star 


Start HJD 


End HJD 


Number of 


Exposure 


Overhead 


Peak 


Telescope/Instrument 




(2450000+) 


(2450000+) 


exposures 


time (s) 


time (s) 


S/N 


(observing mode) 


HD9289 


2920.54506 


2920.62881 


111 


40 


25 


90 


VLT/UVES (600 nm) 


HD 12932 


2921.62234 


2921.70532 


69 


80 


25 


90 


VLT/UVES (600 nm) 


HD 19918 


2921.52607 


2921.60905 


69 


80 


25 


100 


VLT/UVES (600 nm) 


HD24712 


3321.65732 


3321.74421 


92 


50 


22 


300 


VLT/UVES (390+580 nm) 


HD 101065 


3071.67758 


3071.76032 


111 


40 


25 


180 


VLT/UVES (600 nm) 


HD 122970 


3069.70977 


3069.79359 


111 


40 


25 


160 


VLT/UVES (600 nm) 


HD 128898 


3073.80059 


3073.88262 


265 


1.5 


25 


250 


VLT/UVES (600 nm) 


HD 134214 


3070.77571 


3070.85848 


111 


40 


25 


260 


VLT/UVES (600 nm) 


HD 137949 


3071.76312 


3071.84598 


111 


40 


25 


350 


VLT/UVES (600 nm) 


HD 201601 


2871.46470 


2871.56295 


70 


80 


42 


80 


SAO 6-m/NES (425-600 nm) 


HD 201601 


2186.70618 


2186.80296 


64 


90 


43 


250 


Gecko/CFHT(654-6660 nm) 


HD 201601 


2186.82456 


2186.92308 


65 


90 


43 


230 


Gecko/CFHT(662-673 nm) 



as 20 s , thus giving a duty cycle of 70-80% for the majority of 
the targets. The signal-to-noise ratio (S/N) of individual spectra 
is between 90 and 350, as estimated from the dispersion of the 
stellar fluxes in the line-free regions. Detailed description of ob- 
servations for each target is presented in Table [TJ Columns of 
the table give the stellar name, heliocentric Julian dates for the 
beginning and the end of spectroscopic monitoring, number of 
observations, exposure and overhead times, peak signal-to-noise 
ratio of individual spectra, and information about the telescope 
and instrument where the data were obtained. 

The red arm of the UVES spectrometer was configured to 
observe the spectral region 4960-6990 A (central wavelength 
6000 A). The wavelength coverage is complete, except for a 
100 A gap centred at 6000 A. Observations were obtained with 
the high-resolution UVES image sheer (slicer No. 3), provid- 
ing an improved radial velocity stability and giving maximum 
resolving power for A/AA « 1 15 000. 

All spectra were reduced and normalised to the continuum 
level with a routine especially developed by one of us (DL) for a 
fast reduction of spectroscopic time-series observations. A spe- 
cial modification of the Vienna automatic pipeline for echelle 
spectra processing (Tsymbal et al. 2003 ) was developed. All bias 
and flat field images were median-averaged before calibration. 
The scattered light was subtracted by using a 2-D background 
approximation. For cleaning cosmic ray hits, we used an algo- 
rithm that compares the direct and reversed spectral profiles. To 
determine the boundaries of echelle orders, the code used a spe- 
cial template for each order position in each row across disper- 
sion axes. The shift of the row spectra relative to template was 
derived by a cross-correlation technique. Wavelength calibration 
was based on a single ThAr exposure, recorded immediately af- 
ter each stellar time series. Calibration was done by the usual 2- 
D approximation of the dispersion surface. An internal accuracy 
of 30-40 ms~' was achieved by using several hundred ThAr 
lines in all echelle orders. The final step of continuum normali- 
sation and merging of the echelle order was carried out by trans- 
formation of the flat field blaze function to the response function 
in each order. 

The global continuum normalisation was improved by iter- 
atively fitting a smoothing spline function to the high points in 
the average spectrum of each roAp star. With this procedure we 
corrected an underestimate of the continuum level, unavoidable 
in analysis of small spectral regions of the crowded spectra of 
cool Ap stars. Correct determination of the absolute continuum 
level is important for retrieving unbiased amplitudes of radial ve- 



locity variability when the centre-of-gravity method is used. In 
addition to the global continuum correction, spectroscopic time 
series were post-processed to ensure homogeneity in the contin- 
uum normalization of individual spectra. Extracted spectra were 
divided by the mean, the resulting ratio was heavily smoothed, 
and then it was used to correct continua in individual spectra. 
Without this correction a spurious amplitude modulation of pul- 
sation in variable spectral lines may arise due to an inconsistent 
continuum normalisation. 

The red 600 nm UVES dataset was complemented by the 
observations of HD24712 obtained on November 11, 2004 in 
the DDT program 274.D-5011. Ninety-two time-resolved spec- 
tra were acquired with the UVES spectrometer, configured to 
use the 390+580 nm dichroic mode (wavelength coverage 3300- 
4420 and 4790-6750 A). A detailed description of the acquisi- 
tion and reduction of these data is given by Ryabchikova et al. 
(120071) . 

For the roAp star HD 201601 (y Equ), we analysed 70 
spectra obtained on August 19, 2003 with the NES spectro- 
graph attached to the 6-m telescope of the Russian Special 
Astrophysical Observatory. These echelle spectra cover the re- 
gion 4250-6000 A and have typ ical S /N of » 80. The data were 
recorded by Kochukhov et al. (2004), who searched for rapid 
magnetic field variability in y Equ. We refer readers to that paper 
for the details on the acquisition and reduction of the time-series 
observations at the SAO 6-m telescope. 

Post-processing of the echelle spectra of HD24712 and 
HD 201601 was done consistently with the procedure adopted 
for the main dataset. In addition, the time-resolved observa- 
tions of HD 201601 obtained in 2001 using the single-order f/4 
Gecko coude spectrograph with the EEV1 CCD at the 3. 6-m 
Canada-France-Hawaii telescope were used. Observations have 
a resolving power of about 115 000, determined from the widths 
of a number of ThAr comparison lines. The reduction is de- 
scribed in Kochukhov et al. ( 120071) . 

3. Fundamental parameters of programme stars 

Fundamental parameters of the programme stars are given 
in Table [2] For six stars, effective temperatures T e g, surface 
gravity logg, and mean surface magnetic fields (B s ) were 
taken from the literature. For 4 remaining stars, HD9289, 
HD 12932, HD 19918, and HD 134214, atmospheric parameters 
were derived using Stromgren photometric indices (Hauck & 
Mermilliod 1 1 9981) with the calibrations by Moon & Dworetsky 



Table 2. Fundamental parameters of target stars. 



HD 


Other 


Tefi 


l°g£ 


v e sin i 


(B s ) 


P 


Reference 


number 


name 


(K) 




(km s" 1 ) 


(kG) 


(min) 




9289 


BW Cet 


7840 


4.15 


10.5 


2.0 


10.522 


this paper 


12932 


BN Cet 


7620 


4.15 


3.5 


1.7 


11.633 


this paper 


19918 


BTHyi 


8110 


4.34 


3.0 


1.6 


11.052 


this paper 


24712 


DO Eri, HR 1217 


7250 


4.30 


5.6 


3.1 


6.125, 6.282 


Ryabchikova et al. (|1997> 


101065 


Przybylski's star 


6600 


4.20 


2.0 


2.3 


12.171 


Cowley et al. l2000T> 


122970 


PPVir 


6930 


4.10 


4.5 


2.3 


11.187 


Rvabchikova et al. d2000l 


128898 


a Ch- 


7900 


4.20 


12.5 


1.5 


6.802, 7.34 


Kupka et al. (|1996J 


134214 


ill Lib 


7315 


4.45 


2.0 


3.1 


5.690 


this paper 


137949 


33 Lib 


7550 


4.30 


<2.0 


5.0 


8.271,4.136, 9.422 


Ryabchikova et al. (|2004l 


201601 


y Equ 


7700 


4.20 


<1.0 


4.1 


12.20 


Rvabchikova et al. J2002l 



d!985t and by Napiwotzki et al. d 19931 1 as implemented in the 
TEMPLOGG code (Rogers 1995 ). In addition Geneva photomet- 
ric indices (Burki et al. 2005 Q with the calibration of Kiinzli 
et al. ( 1 19971 1 was used for effective temperature determination. 
The colour excesses were estimated from the reddening maps by 
Lucke (119781 1. In Table [2] we present average values of the ef- 
fective temperatures derived with three different calibrations. A 
typical dispersion is +150 K. 

In all stars but HD 101065, rotational velocities were es- 
timated by fitting line profiles of the magnetically insensitive 
Fei 5434.5, 5576.1 A lines. Magnetic spectral synthesis code 
SYNTHMAG (Piskunov [T999l Kochukhov l2006al l was used in our 
calculations. Atomic parameters of spectral lines were extracted 
from the VALD (Kupka et al. 119991 1 and DREAM (Biemont et al. 
1999) databases, supplemented with the new oscillator strengths 
for Lan (Lawler et al. HEP, N dn (Den Hartog et al. [2003), 
Ndm (Ryabchikova et al. 20063, Smu (Lawler et al. 120061) . and 
Gdn (Den Hartog et al. 2006). We confirmed rotational veloci- 
ties derived previously for HD 24712, HD 122970, HD 128898, 
HD 137949, and HD 20 1601. High spectral resolution of the 
present data allows us to improve the value of the projected rota- 
tional velocity for HD 101065 (Przybylski's star) using partially 
resolved Zeeman patterns in numerous lines of the rare-earth ele- 
ments. The value of the magnetic field modulus, 2.3 kG (Cowley 
et al. 120001 ). was confirmed. 

We note that the derived values of v e sin i in very sharp-lined 
roAp stars depend strongly on the Fe stratification, which is 
known for some of them (e.g. y Equ, see Ryabchikova et al. 
120021 ) and is typically characterised by the Fe overabundance 
below logT5ooo = -1 and Fe depletion in the outer layers. 
Rotational velocities in Table [2] are derived for a homogeneous 
Fe distribution. However, in all programme stars we see observa- 
tional evidence of Fe stratification, such as anomalous strength 
of the high-excitation lines compared to the low-excitation lines 
(see Ryabchikova et al. 2003 ). We tested one of the most peculiar 
stars, 33 Lib, for possible influence of stratification on the de- 
rived v e sin i and found that including stratification may decrease 
the inferred rotational velocity to 1.5 kms -1 in comparison to 
2.5-3.0 km s _1 obtained with a homogeneous Fe abundance dis- 
tribution. The lower value of v e sin i seems to fit the observed 
resolved and partially resolved Zeeman components better. 

In three programme stars, HD9289, HD 12932, and 
HD 19918, mean magnetic modulus was estimated for the first 
time from differential magnetic broadening/intensification. Two 
spectral regions with a pair of lines having different magnetic 
sensitivity were synthesised for a set of magnetic field strengths. 
The first region contain well-known Fen 6147.7 A (g e ff=0.83) 

1 http : //obswww . unige . ch/gcpd/phl3 . html 



and 6149.3 A (g eff =1.35) lines, while Fei 6335.3 A (g eff =1.16) 
and 6336.8 A (g e ff=2.00) lines were analysed in the second spec- 
tral region. The derived values of the magnetic fields were con- 
firmed by fitting other magnetically sensitive lines, for example, 
Fen 6432.7 A (g eff = 1.82) and Eun 6437.6 A (g eS = 1.76). 



4. Radial velocity measurements 

To perform a meaningful study of the pulsational amplitudes 
in spectral lines of different chemical elements/ions, one has 
to be very careful in the choice of lines for pulsation mea- 
surements. For this purpose we have synthesised the observed 
spectral region for each star with the model atmosphere param- 
eters and magnetic field values from Table [2] Abundances for 
HD 24712, HD 101065, HD 122970, HD 128898, HD 137949, 
and HD 201601 were taken from the papers cited in the last col- 
umn of Table|2l For the remaining four stars, a preliminary abun- 
dance estimate was obtained in this paper. 

The radial velocities were measured with a centre-of-gravity 
technique. We used only unblended or minimally blended lines. 
In some cases where the line of interest was partially overlap- 
ping with the nearby lines, only the unblended central part of the 
line was considered; therefore, some lines were not measured 
between the continuum points. This usually leads to lower pul- 
sation amplitudes if we have strong variations in the pulsation 
signal across a spectral line (see y Equ - Sachkov et al. 2004 
and HD 99563 - Elkin et al. 2005 ). Bisector radial velocity mea- 
surements were performed for Ha core and for a subset of Y n, 
Eun, Ndn, Ndm, Pru, Prin, Tbm, and Thra spectral lines. The 
pulsational RV variability of the two strongest Th m lines at AA 
5376.13 and 6599.48 A was investigated for several stars in our 
sample for the first time. 

A detailed frequency analysis of the RV data for 8 stars from 
our sample was carried out by Kurtz et al. (120061) . These au- 
thors used the same observations as we do and, despite the rela- 
tively short time span (2 hours) of the spectroscopic time series, 
they claimed to resolve several frequencies for each star and to 
find amplitude modulation that was not observed in photometry. 
Although a detailed frequency analysis is not the primary goal of 
our paper, we did repeat time-series analysis for all lines. First 
we applied the standard discrete Fourier transformation (DFT) 
to the RV data. The period corresponding to the highest pulsa- 
tion amplitude value was then improved by the sine-wave least- 
square fitting of the RV data with pulsation period, amplitude, 
and phase treated as free parameters. This fit was removed from 
the data and then Fourier analysis was applied to the residuals. 
This procedure was repeated for all frequencies with the S/N 
above 5. 



To verify our analysis against the results of Kurtz et al. 
(2006 ), we applied it to the set of Prin lines (to each line sep- 
arately and to the average Prm RV) and compared in detail 
the resulting solution for HD 134214. Our frequency solution 
agrees perfectly with Kurtz et al. (2006), although our RV am- 
plitudes are systematically lower. We believe that the primary 
reason for this discrepancy, which is present for other stars as 
well, is our systematically higher continuum placement, leading 
to lower pulsation amplitude when the centre-of-gravity method 
is used. As explained in Sect. [2] we have rectified observations 
with a spline-fit over wide wavelength regions, whereas Kurtz et 
al. (f2006 ) have probably assigned continuum points to the high 
spectrum points in the immediate vicinity of the line considered. 

For each spectral line, we estimated the probability that the 
detected periodicity is not due to noise (Home & Baliunas 1986). 
Then, for the lines with probabilities higher than 0.999, we cal- 
culated a weighted average value of the pulsation period. All 
amplitudes and phases were then recalculated keeping the pe- 
riod fixed. This information is summarised in Table [3] (available 
online only) for all lines studied in each star. Table [3] contains 
amplitudes (first line) and phases (second line), together with the 
corresponding errors for the dominant pulsation period. When 
no pulsation signal was detected, we give the formal amplitude 
solution for the fixed period without phase information. If a star 
had one dominant period, the RV analysis was done with this 
period. In the cases of more than one dominant period, a simul- 
taneous fit was performed with up to three periods. RV variation 
was approximated with the expression 

3 

(V) = V (t - 1 ) + Yj y i cos < 2?r [( f - f o)/^ - fi]l (1) 
i=l 

Here the first term takes possible drift of the spectrograph's zero 
point into account. For all stars but HD 24712, HJD of the first 
exposure of the star at a given night was chosen as a reference 
time t . For HD 24712, HJD=2453320.0 was used as a reference 
time. Both V, and <fi are, respectively, amplitude and phase of 
the RV variability with the ith period (/ max = 3). With the minus 
sign in front of <fj, a larger phase corresponds to a later time of 
RV maximum. This phase agreement is natural when discussing 
effects of the outward propagation of pulsation waves in the at- 
mospheres of roAp stars. 

The periods employed for determining of the RV amplitudes 
and phases are given in the seventh column of Table [2] When 
more than one period was inferred in the fitting procedure, the 
main period for which pulsation amplitudes and phases are given 
in Table [3] (online material) is highlighted with italics. 

5. Bisector measurements 

Radial velocity analysis was complemented by studying bisec- 
tor variability. Figure Q] shows bisector amplitudes (left panel) 
and phases (right panel) across the core of Ha line calculated 
with the pulsation periods from Table [2] Where two periods or 
the main period and its first harmonic are resolved, simultane- 
ous fit with two frequencies was done. In all programme stars, 
an increase in bisector amplitude by two or more times from the 
transition region to the deepest part of the core is observed. Thus 
we supported the conclusion made by Kurtz et al. ((2005a). This 
change is gradual in all but two stars: HD 9289 and HD 19918, 
which have the lowest S/N. It is, therefore, unclear if small jumps 
of RV amplitude across the core are real or caused by the low S/N 
of the spectroscopic data. RV changes are accompanied by tiny 
phase changes. Only in four stars do phase changes exceed the 



error bars. These are the stars with the shortest pulsation periods 
close to the acoustic cut-off frequencies: HD 24712, HD 128898, 
HD 132214, and HD 137949 (33 Lib). The NLTE calculations 
show that, in the atmosphere of the star with T s s between 7000 
and 8000 K, the Ha core is formed at -5 5; logrsooo ~ -2 
(Mashonkina, private communication). It is applied to a normal 
atmosphere; however, a core-wing anomaly is present in all pro- 
gramme stars (Cowley et al. 2000 ), which was attributed to a pe- 
culiar atmospheric stratification containing a region of increased 
temperature below log T5000 = -4 (Kochukhov et al. 2002). This 
change in the atmospheric structure may lead to an upward shift 
in formation depth of the base of the Ha core (Mashonkina, pri- 
vate communication). 

More care should be taken in the choice of metal spectral 
lines for the bisector measurements than for the centre-of-gravity 
RV analysis. Blends with non-pulsating lines change the run 
of bisector pulsation amplitude across the line profiles, while 
blends with pulsating lines (weak REE lines in far wings, for ex- 
ample) may change both amplitudes and phases. These artifacts 
often lead to the wrong conclusion about roAp pulsational be- 
haviour, especially when these conclusions are based on analysis 
of only one line. For instance, the well-known Nd m 6145.068 A 
line, analysed in many pulsation studies, is blended with the 
Sii 6145.016 A line and with the Cen 6144.853 A. The latter 
feature makes a non-negligible contribution to the spectrum of 
Przybylski's star. The strongest Ndm 5294.11 A line is blended 
with the Fei 5293.96 A line in the blue wing and with the high- 
excitation Mnn 5294.32 A line in the red wing. Both lines are 
normally weak, but they are strengthened in the stratified atmo- 
spheres of roAp stars with the effective temperatures higher than 
7000-7200 K. This blending results in a drop of the bisector ve- 
locity amplitudes starting from some intensity points in the line 
profile. Zeeman splitting may also change velocity amplitude 
distribution across the line profile. Figure |2] illustrates the influ- 
ence of blends and Zeeman effect on the bisector amplitudes and 
phases of Ndn- m and Prm lines in HD 134214. Pairs of Ndm 
lines with similar intensities and Zeeman structure have identi- 
cal bisector phases independent of any blending. Thus, phases 
may be considered as a more reliable indicator of the pulsational 
characteristics. One may see that, up to the residual intensity 
0.85, the Nd m 5294 A line shows the same bisector velocity am- 
plitudes as another strong Nd m line at 5204 A. However, after 
this intensity point, the two amplitude curves start to deviate: 
in the Nd m 5294 A line, the amplitude drops towards the line 
wings, probably due to blends. Also blends decrease bisector 
amplitudes in the Ndm 6 145 A compared to the Ndm 6327 A 
line, which has the same intensity and Zeeman structure. Three 
out of the four Prm lines shown in Fig.|2]are practically free of 
blends in most of the programme stars, and additionally, these 
lines have identical Zeeman splitting, which results in similar 
bisector pulsational behaviour. 

As shown in Section|6] Tb m lines exhibit a particularly in- 
teresting pulsation behaviour, so their blending should be in- 
vestigated in detail. The Tbm lines are not as strong and nu- 
merous as Prm and Ndm. From 10 lines only the strongest 
one, Tbm 5505 and three others, 5847, 6323, and 6687 A, 
may be used for bisector measurements. Most often we mea- 
sured the Tbm 5505.408 A line. This line is blended by the 
Cen 5505.204 A in the blue wing and by Cen 5505.588 A in 
the red wing. The Tbm A 6323.619 A line is free of blends but 
may be partially blended with the weak atmospheric O2 line at 
A 6323.75 A. Figure [3] shows a spectral region around the Tb m 
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Fig. 1. Bisector measurements in the Ha core for the programme stars. The RV amplitudes are shown in the left panel and pulsation 
phases (based on the main periods from Table 13 in the right panel. 
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Fig. 2. Bisector measurements as a function of the Nd (left panels) and Pr (right panels) lines' residual intensity in HD 134214. Top 
panels represent RV amplitudes, and bottom panels show pulsation phases calculated with the main period from Table [2] 



5505.408 A line in HD 101065, where the blending problem is 
the most severe. For comparison we show the region with the 
Tbm 6323.619 A line; therefore, wavelength scale is given in 
km s _1 relative to the centre of each line. We also show synthetic 
spectrum calculated with the highest possible Ce abundance and 
magnetic field strength appropriate for HD 101065 (see Table©. 
Although the full intensity of the calculted feature is equal to the 



intensity of the observed one, for demonstration of the blend- 
ing effects we did not convolve synthetic spectrum with the in- 
strumental and rotational profiles. Note that the red wing of the 
Tbm 6323.619 A line is blended with the atmospheric O2 line 
6323.75, therefore no bisector measurements were carried out 
above intensity level 0.9. The blue wing asymmetry in both Tb m 
lines is rather an effect of the hyperfine splitting, than the real 




Wavelength, km/s Residual intensity Residual intensity 

Fig. 3. Tb in lines in HD 101065. Left panel represents the observed and calculated line profiles, and the two next panels show the 
RV amplitudes and pulsation phases determined for the intensity levels indicated in the corresponding line profiles by the horizontal 
lines. 



influence of blends, because the same asymmetry is observed in 
other stars where the Ce n contribution is negligible. 

Taking into account the similarity of the RV amplitudes and, 
in particular, phases in the bisector measurements of both lines 
we conclude that even in HD 101065 the blending does not affect 
seriously the measured pulsation characteristics. 

6. Phase-amplitude diagrams 

Empirical understanding and detailed theoretical interpretation 
of the pulsation phenomenon in roAp stars requires construction 
of the vertical cross-section of pulsation modes. Chemical strat- 
ification enables this difficult task by separating formation lay- 
ers of the spectral lines of different chemical elements far apart, 
thus allowing us to resolve the vertical dimension of pulsating 
Ap-star atmospheres, as can be done for no other type of pulsat- 
ing stars. However, the key information needed for such vertical 
pulsation tomography - formation depth of REE spectral lines - 
is difficult to obtain. Early studies (e.g., Kanaan & Hatzes 1998) 
have suggested that the line intensity may be used as a proxy 
of the relative formation heights. However, it is now understood 
that the chemical stratification effects are dominant in the atmo- 
spheres of cool Ap stars and, therefore, formation region of weak 
lines of one element is not necessarily located deeper than the 
layer from where the strong lines of another element originate. 
This is why a comparison of the intensities of different pulsating 
lines is only meaningful for the absorption features of the same 
element/ion. 

The study by Ryabchikova et al. (2002) presented the first 
comprehensive analysis of the vertical stratification of light, 
iron-peak, and rare-earth elements in a roAp star, demonstrat- 
ing that only by taking height-dependent abundance profiles into 
account can one calculate correct formation depths of pulsat- 
ing spectral lines and meaningfully interpret the results of time- 
resolved spectroscopy. In a later paper, Mkrtichian et al. (2003 ) 
study the vertical profiles for p-mode oscillations in 33 Lib 
(HD 137949) using pulsation centre-of-gravity measurements 
for individual lines. To establish the vertical atmospheric coor- 
dinate, the authors used the Wa - log t transformation scale in 
33 Lib's atmosphere by assuming homogeneous elemental dis- 
tribution. This simplified approach was sufficient for detecting 
the node in the Ndn, Ndm line-forming region - a discovery 
later confirmed by Kurtz et al. d2005bl l and by the present paper. 
However, as Mkrtichian et al. (120031 ) themselves acknowledge, 



their method is too simplified because it ignores the vertical 
and horizontal inhomogeneous distribution of chemical abun- 
dances. Indeed, their technique cannot be used for simultaneous 
interpretation of the pulsational variability of different elements 
since Fe, as well as Ca, Cr and all REE elements, are strongly 
stratified in the atmosphere of 33 Lib and other roAp stars (see 
Ryabchikova et al. 2004), and stratification is not the same for 
different elements. Mkrtichian et al. (120031 ) suggest that in future 
studies some of these difficulties can be alleviated by the "multi- 
frequency tomographic approach" that will use the phase and 
amplitude profiles of different oscillations modes to establish a 
link to the atmospheric geometric depth scale. 

In general, the detailed pulsation tomography analysis of 
roAp stars should be based on sophisticated atmospheric mod- 
elling including chemical stratification, NLTE, magnetic field, 
and eventually, pulsation effects on the shape and intensity of 
spectral lines with substantial RV variability. Such modelling is 
very demanding in terms of the quality of observations, required 
input data, and computer resources. This is why only two roAp 
stars, y Equ and HD 24712, have been studied with this method 
up to now (Mashonkina et al. 2005; Ryabchikova et al. [2006b). 
Here we suggest a different approach to the pulsation tomog- 
raphy problem. In the framework of the outward propagating 
magnetoacoustic wave, one expects a continuous amplitude ver- 
sus phase relation for pulsation modes. Lines showing later RV 
maximum should originate higher in the atmosphere. Thus, new 
insight into the roAp pulsation modes structure can be obtained 
by inferring and interpreting the trend in pulsation velocity am- 
plitude as a function of pulsation phase. Such phase-amplitude 
diagrams offer a possibility of tracing the vertical variation in the 
mode structure without tedious assignment of the physical depth 
to each pulsation measurement. 

We note that this paper is not the first one to use the phase- 
amplitude diagrams for roAp stars. Baldry et al. (1998]) and 
Baldry & Bedding (2000) produced similar diagrams in the pul- 
sational studies of a Cir and HR 3831. However, their diagrams 
described pulsational properties not for individual lines with a 
proper identification but for small spectral regions often con- 
taining several spectral lines, sometimes without correct iden- 
tification. As a result, the most important information on wave 
propagation was lost (see Sect. 16.51 ). 

For each programme star we have produced phase-amplitude 
diagrams using a set of lines of representative chemical species, 
including Y n, Eu n, the Ha core, La n, Dy n, Dy m, Nd n, Nd m, 
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Fig. 4. Amplitude-phase diagrams for the pulsation centre-of-gravity measurements (left panel) and bisector measurements (right 
panel) in HD 101065. 



Pr ii, Pr m, Tb m, and Th m. We considered velocity amplitudes 
and phases derived with both the centre-of-gravity RV measure- 
ments and the bisector analysis. The line blending varies signif- 
icantly from star to star due to different effective temperatures 
and chemical anomalies, therefore we cannot use identical set 
of spectral lines for all stars. For each star we tried to employ 
unblended or minimally blended lines, in particular, for the bi- 
sector measurements. In addition, the line broadening expressed 
in the terms of macroturbulent velocity was estimated for a few 
representative lines of each chemical species. This information 
is useful for assessing the isotropic velocity component at the 
formation heights of REE lines. 

Below we present the results for individual stars. 

6.1. HD 101065 (Przybylski's star) 

The atmosphere of HD 101065 is known to be very rich in REE 
and underabundant in most other elements, including the Fe- 
peak species. Due to the low effective temperature, slow rota- 
tion, and abundance anomalies, most lines in the spectrum of 
HD 101065 are strong and sharp. As a result, one can achieve 
impressive accuracy of pulsation measurements, typically ~5- 
8 ms~', for moderately strong lines. Although it is not easy to 
find unblended Fe lines in the forest of the REE lines, we man- 
aged to measure a few of them and to detect no pulsational vari- 
ability above 30 ms~'. Pulsation amplitudes at the level of 6 to 
17 m s _1 were detected in the cores of very strong Ba n AA 6142, 
6496 A lines, which gives us an idea about pulsation amplitudes 
close to the photospheric layers. 

The centre-of-gravity and bisector phase-amplitude dia- 
grams for HD 101065 are shown in Fig. [4] The centre-of-gravity 
measurements give us a general idea of the pulsation wave prop- 
agation in roAp atmosphere. Pulsations are characterised by the 
amplitude ~40 ms~' at the levels of the Yn, La a, and strong 
Eu ii line formation. Then they pass, with the gradually increas- 
ing amplitude, the layers where the Dy n, Dy m, Nd n, Pr n, Nd m, 
the Ha core and Prm lines are formed. The phase does not 
change by more than 0.15 of the pulsation period (less than 1 
radian) from Eun to Prm lines. After that, a rapid change in 
pulsation phase occurs in Tbin and, next, in Thin lines. We 
measured two unblended Thm lines at A 5376.13 and 6599.48 A. 

HD 101065 is the only star in our sample that shows similar 
pulsation signatures (bisector amplitude versus pulsation phase) 
for the Ha core, Ndn, Ndm, Prn and weak Prm lines. For 



all these features we find an increase in both pulsation ampli- 
tude and phase from the line wings to the line core. Kurtz et 
al. (2005a) have obtained similar amplitudes and phases for the 
Ha core and, based on the bisector measurements of the single 
Nd in 6 145 A line, they argue that the formation layers of Nd m 
lines start above the layers of formation of the deepest part of 
the Ha core. Our Fig. [4] provides strong evidence that most of 
REE lines, including Pr n and weaker Pr m lines, are formed in 
the same layers as the Ha core. 

In addition we analysed the broadening of pulsating lines. 
All low-amplitude Ce n lines do not exhibit any extra broaden- 
ing above the adopted value of the projected rotation velocity. 
All Eu n lines and Dy and Nd lines in both ionisation states show 
an extra broadening corresponding to V ma cro between 2.5 and 
4 km s _1 . Macroturbulent velocity grows from 5 km s in weak 
Pr m 5765 A line to 7-8 km s _1 in the strong Pr m lines AA 5300, 
6707 A. We need V macro = 10 km s to fit the observed profiles of 
the unblended Tb in 6832 A and Thm 5376 A lines. It is worth 
noting that this rapid increase in the line broadening accompa- 
nies the change in the line profile variability pattern from the 
normal symmetric shape to the blue-to-red running waves, which 
become noticeable for the strongest Pr in lines and, in particular, 
for Tb m lines (Kochukhov et al. 2007). An extra broadening of 
the Pr m lines cannot be attributed to hyperfine structure. 

Considering the amplitude-phase diagrams, we see a change 
in the pulsation behaviour for Tb in lines and, even stronger, for 
the Thin 5376.13 A line. In contrast to the variability of the 
Ha core and Nd lines, we observe a rapid increase in the am- 
plitude and in the phase from the line core to line wings in the 
lines of Tbm and Thm (see Fig. [T2]- Online only). This phe- 
nomenon cannot be explained by the pulsation wave propagating 
outwards. 

In general, over the large fraction of the Przybylski's star 
atmosphere, up to the layers of the Ndm and weak Prm line 
formation, pulsations are represented well by non-radial non- 
adiabatic /9-modes (Saio 2005) with an amplitude increasing out- 
ward. Above these layers, pulsation characteristics change in a 
way that is not forseen by the theoretical models. 

6.2. HD 122970 

Pulsational observations of this star were obtained at the ro- 
tational phase 0.75, between the maximum and minimum of 
the magnetic field (see Fig. 3 in Ryabchikova et al. 2005a). 
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Fig. 5. Bisector measurements as a function of lines residual in- 
tensity in HD 122970. The left panel represents RV amplitudes, 
right panel shows pulsation phases calculated with the main pe- 
riod from Table [2] 



The centre-of-gravity and bisector phase-amplitude diagrams for 
HD 122970 are presented in Fig. [13] (Online only). For a few 
chosen lines we also show the variability of the pulsational char- 
acteristics along the line profiles in Fig. [5] 

The overall pulsational behaviour in the atmosphere of this 
star is similar to HD 101065: low-amplitude pulsations are seen 
in Yn and Eun lines, the amplitude increases slightly in the 
layers where the Ha core forms, and then it increases further 
in the Dy and Nd line formation layers. A rapid growth of the 
amplitude occurs in the region of the Pr line formation. Up to 
these layers, pulsations have a standing-wave character, with al- 
most constant phases for all lines of Y n, Eu n, La n, Dy, Nd, 
and Pr. The phase shifts become noticeable for Tb m and Th in 
lines, as in HD 101065, and, again, in these lines pulsation am- 
plitude and phase increase from the line core to line wing, con- 
trary to predicted effect of the outward-running wave. While in 
HD 101065 the regions of the Ho- core, Nd, and Pr line forma- 
tion are the same or, at least, partially overlap, in HD 122970 the 
phase-amplitude diagrams show a significantly more stratified 
atmospheric structure. 

Ryabchikova et al. (2000) noted that the projected rotational 
velocities derived from the REE lines are lower than those ob- 
tained from Fe-peak lines. More precise present observations 
confirm this result. We determined v e sin i » 3.5 km s 'from the 
REE line fitting, whereas the Fe-peak elements show broaden- 
ing corresponding to v e sin i = 4.5 km s" 1 . A 3.877 d rotation pe- 
riod obtained by Ryabchikova et al. (2005a) for HD 122970, to- 
gether with the low v e sin i, suggests a nearly pole-on geometry. 
Therefore, the difference in the apparent rotational velocities of 
the Fe-peak and REE lines may indicate a concentration of the 
REEs near the visible stellar rotation pole. No extra broaden- 
ing is required to fit the Nd n, Nd m, Pr n, and weak Pr m lines, 
while V m acio = 2-3 kms -1 and 6 kms~' are needed to fit the 
strong Prm lines and the Thin 5376 A line, respectively. The 
Yn lines, which show the smallest detectable RV amplitudes, 
have the same pulsation phases and extra broadening as strong 
Prm lines. 

The lines of the elements with Z < 38 do not show any pulsa- 
tion amplitude above 7-20 m s _1 . Two measured Zr n lines have 
amplitudes of « 200 m s _1 . 



6.3. HD 24712 

According to the ephemeris given by Ryabchikova et al. (2005a), 
our pulsational observations were obtained at the rotational 
phase 0.944, i.e. close to the magnetic maximum. Detailed ob- 
servational analysis of the spectroscopic pulsation signatures in 
HD 24712 was recently presented by Ryabchikova et al. ( 120071 ). 
HD 24712 is the only roAp star for which NLTE line formation 
calculations of Ha, Pr, and Nd were carried out (Mashonkina et 
al. 2005; Ryabchikova et al. 2006b). The observed intensities of 
the Nd and Pr lines in the first and the second ionisation states are 
explained by the stratified atmosphere with a step-like enhance- 
ment of these elements just above the NLTE formation depth of 
the Ho- core (logT 5 ooo = -4). Up to this atmospheric level, the 
observed distribution of pulsation phases agrees rather well with 
the predictions of the non-adiabatic pulsation models (see Fig. 1 
in Kochukhov 2006b and Fig. 3 in Sachkov et al. 120061) . which 
is supported by the amplitude -phase diagrams for HD 24712 
shown in Fig. [14] (Online only). In contrast to the standing-wave 
dominated behaviour of HD 101065 and HD 122970, we see a 
running wave in HD 24712 from the low atmosphere, where the 
first pulsation signal is detected. Due to the availability of the 
simultaneous spectroscopic and photometric monitoring, the se- 
quence of the pulsation phase changes is established with high 
accuracy in this star (Ryabchikova et al. 2007). Pulsation follows 
the same order as in HD 101065 and in HD 122970: Eu, La, Ha 
core, Nd, Pr, Tb lines. 

The phase distribution with the optical depth, as well as our 
amplitude-phase diagrams, shows the phase shifts between Ha 
and Nd lines and between Pr and Nd lines. The NLTE calcula- 
tions do not predict a significant difference in the location of the 
Pr and Nd line formation regions, therefore the observed phase 
shift is probably caused not by the vertical distribution of chem- 
ical elements, but by the horizontal abundance inhomogeneities. 
According to the pulsation phases, Nd lines have to be formed at 
the same layers as the deeper parts of the Ha core. At the same 
time, we caution that the existing NLTE calculations of Ha, Pr, 
and Nd lines are very preliminary, and they do not include poten- 
tially important effects of stratified abundance distribution on the 
atmospheric structure and neglect an influence of the magnetic 
desaturation on the line formation. 

The bisector velocity amplitudes of HD 247 12 are either con- 
stant along the line profile, or they slightly increase from the 
line wings to the line core, except for the Tb in and, possibly, 
Yn lines, in which bisector amplitude increases from the line 
core to the line wings (Fig. [15]- Online only). Pulsation phases 
are nearly constant along the profiles of most lines, but show a 
tendency to grow from the line centre to the line wings for the 
strongest Prm lines and for Tbm lines. Both Thin lines show 
pulsational amplitudes close to the detection limit. A very in- 
teresting result is obtained for Yn lines. Yttrium is the only 
element with Z<40 whose lines have a measurable pulsation 
amplitude in HD 24712 (see Ryabchikova et al. 120071 ). These 
lines are strong enough to provide precise centre-of-gravity and 
bisector measurements, therefore pulsation phases can be de- 
termined with high accuracy. While RV amplitudes do not ex- 
ceed 50-100 ms~' for the strongest Yn 5087 A line, the bisec- 
tor phase distribution across the line profile coincides with the 
phase distribution across the profile of the strongest Pr m 5300 A 
line, which has an RV amplitude above 300 m s _1 . The estimated 
depth of formation of Y n lines in chemically homogeneous at- 
mosphere lies around log T5000 = -1 -5 to -2.5, where the core of 
Ha line starts forming. Note that Doppler imaging of HD 247 12 
shows a similarity between yttrium and REE surface distribu- 



tions (Liiftinger et al. 120061 1. The observed phase difference be- 
tween Y ii and the Ha and the coincidence of the phases in Y n 
and strong Pr m lines, on the one hand, and the large difference 
in amplitudes, on the other, is difficult to explain. We will see 
that the same pulsational behaviour of Y n lines is observed for 
some other roAp stars. 

Although the projected rotational velocity of HD 24712 is 
rather high, spectral synthesis clearly shows that we need an 
extra broadening, equivalent to Vmacro = 6 kms -1 , to fit strong 
Prm lines, as well as Tbm and Thin lines, whereas no extra 
broadening is required for Nd lines. Interestingly, one needs 
Vmacro = 4 km s -1 to reproduce the profile of the Y n 5087 A line. 

6.4. HD 134214 

This star has one of the shortest pulsation periods, 5.7 min, 
which is close to the acoustic cut-off frequency calculated for 
HD 134214 by Audard et al. (IT9981 . The centre-of-gravity and 
bisector amplitude-phase diagrams are shown for HD 134214 in 
Fig. Q2] (Online only). The amplitude-phase diagrams for this 
star are similar to those for HD 24712. In both stars we ob- 
serve the running wave from the regions of Eu n line formation 
to Tb in, Th m line-forming regions. In HD 1 342 14, all lines pul- 
sating with high amplitude show an increase both in amplitude 
and in phase from the line core to the line wings (see Fig.|2|. 

The low rotational velocity allows us to partially resolve 
components of the Zeeman split lines, which makes the line pro- 
file analysis easier and, in particular, allows us to estimate an 
extra broadening accurately. An extra broadening of ~1 kms -1 
is needed to fit the low-amplitude pulsating lines of Eu n. Then 
the inferred V macl0 increases gradually: 4 kms -1 (Ndn, Dym), 
5-6 kms -1 (Ndm), 8-10 km s" 1 (Prm), 10 kms -1 (Tbm), and 
more than lOkms -1 (Thin). Just as in HD 24712, Yn lines have 
small detectable RV amplitudes, but the same bisector phases as 
strong Prm lines do. We measured V ma cro ~6 kms -1 for the Yn 
5087 A line. 

In both stars the bisector amplitude-phase diagrams for Nd n 
and weaker Ndm lines are overlapping, which supports the 
model of the Nd stratification in roAp atmospheres proposed by 
Mashonkina et al. (2005). According to this model the strongest 
lines of Nd n are formed at the same layers as weak Nd m lines, 
therefore their pulsational characteristics should be similar. 

No definite pulsation signatures are detected in the lines with 
Z < 38 and also in Zr n and Ba n lines. The upper RV amplitude 
limit is 10 ms -1 . 



6.5. HD 128898 (a Cir) 

The main photometric period of a Cir, one of the short-period 
roAp stars, is 6.83 min, which is close to the acoustic cut- 
off frequency. The amplitude-phase diagram for HD 128898 
(Pig. |6j appears to be similar to those constructed for HD 24712 
and HD 134214 and indicates a running wave. Based on an 
amplitude-phase diagram made for short spectral regions con- 
taining groups of unresolved lines, Baldry et al. (1998) sug- 
gested a high-overtone standing wave with a velocity node in 
the atmosphere of a Cir. This contradicts our results obtained 
from analysis of pulsations in isolated lines. Again, similar to 
HD 24712 and HD 134214, Yn lines have the smallest RV am- 
plitudes, but the same phases as the strongest Ndm and Prm 
lines. 

The bisector analysis shows that bisector amplitude rapidly 
decreases from the line centre to the line wings with small 
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Fig. 6. Amplitude-phase diagrams for the pulsation centre-of- 
gravity measurements in a Cir. 



changes in phase. This behaviour may be explained by a com- 
bination of non-radial pulsations and surface chemical inhomo- 
geneity, which clearly manifests itself in the rotational modula- 
tion of the line profiles in a Cir (Kochukhov & Ryabchikova 
2001b). The average spectrum used in the present study cor- 
responds to the rotation phase 0.44 (minimum of the photo- 
metric pulsational amplitude) if we apply the ephemeris from 
Kurtz et al. (1994). All REE lines in our a Cir spectrum ex- 
hibit two-component profiles, which represents another proof of 
chemical spots. 

6.6. HD 12932 



Kurtz et al. (|2006) suggest that HD 12932 is an example of the 
star with a standing wave behaviour in the layers of the Ha 
core formation based on the constancy of the bisector pulsa- 
tion phases. Our measurements are presented in Fig. [TTl (Online 
only). While for Nd n and the Ha core the bisector phases are al- 
most constant for a given line, the phases differ from line to line. 
Phase shifts are small, but reliable, taking the high accuracy of 
the phase determination into account (Fig. [18]- Online only). The 
bisector phase measurements show that, if we attribute a given 
phase to the particular atmospheric height, many Nd lines should 
originate in extremely thin separate layers, which is impossible 
to explain. 

Due to the slow rotation and a rather weak magnetic field, 
HD 12932 is useful for demonstrating the effect of differen- 
tial extra broadening of the pulsating lines. Figure [7] presents a 
comparison between the line profiles of different elements/ions, 
scaled to the same central line depth. The wavelength scale is 
given in kms -1 , with line centre at zero velocity. The two iron 
lines with zero and large Zeeman broadening are shown to as- 
sess the expected magnetic field broadening effect. These lines, 
Fei 5434 A and Fe n 6432 A, do not require any additional broad- 
ening other than Zeeman, instrumental, and rotational ones to fit 
the observed line profiles. The Eu n line broadening is explained 
by the combined Zeeman and hyperfine splitting. Our spectral 
synthesis shows a growth in the Vmacro from Ndn lines (3-5 
kms -1 ), to Ndm lines (6-9 kms -1 ), Prm lines (7-10 kms -1 ), 
and finally, to Tbm lines (10-11 kms -1 ). The stronger line of 
the same ion, the larger Vmacro is required to fit the line profile. 
The Thm lines in HD 12932 are too weak for reliable measure- 
ments. The Y ii 5087 A line exhibits an extra broadening similar 
to Nd ii lines, although no definite pulsation signal was detected 
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all these stars we find first an increase in the RV amplitudes with 
approximately constant phases (standing wave), and after that 
the pulsation wave transforms into a running one. 

As in HD 12932, differential extra broadening is required to 
fit line profiles: V macro = km s (Fe), V macro = 2-3 km s _1 (Y n), 
V ma cro = 4 km s" 1 (Eu n), V macro = 7-8 km s" 1 (Nd n), V macro = 8- 
10 km s _1 (Nd m and Pr ra), and V macl0 = 1 1 km s _1 (Tb ra). Note 
that a smaller broadening is necessary to fit the line core than the 
line wings. Most line profiles show a triangular shape. 

Kurtz et al. (2005a) report a detection of weak RV oscilla- 
tions in the Fei 5434 A line. Our measurements confirm this re- 
sult and also provide the following pulsation RV amplitudes for 
light and iron-peak elements: 30-50 ms~' for Na, Ca, Cr, Fe 
lines; *100ms 1 for Ban lines; 150-200 ms~' for Tin lines. 



Fig. 7. A comparison of the line profiles of different elements in 
HD 12932. Line profiles are scaled to the same central depth. 



in all measured Y n lines, as well as in the lines of Na, Mg, Si, 
Ca, Sc, Cr, Fe, Co, Ni, and Sr. No variability is detected in the 
Fei 5434 A line (RV amplitude 7+15 ms _I ). Thus, we do not 
confirm results by Kurtz et al. d2005a). who claim a 3cr detec- 
tion of pulsations in this line. We looked at 47 other lines of Fei 
and Fen, and we detected a pulsation signal exceeding the 3<x 
significance level only in seven features. For most of these lines 
weak variability is explained by blending with weak REE lines. 
Three out of six measured Ti n lines show a pulsation signal with 
the RV amplitudes 80-100 m s -1 , while two Ti lines did not re- 
veal any variability. All three Ba n lines show RV amplitudes of 
^200 ms- 1 . 



6.7. HD 201601 (y Equ) 

The amplitude-phase diagrams for the centre-of-gravity mea- 
surements in y Equ are shown in Fig. [19] (Online only; left 
panel). Rather low spectral resolution of the y Equ observations 
used for this type of pulsation measurements does not allow 
us to make a precise bisector study. For the bisector analysis 
we used the time-series observations of y Equ carried out with 
the Gecko spectrograph at CFHT (see Kochukhov et al. [2007). 
The bisector amplitude-phase diagrams are displayed in Fig. IT9l 
(right panel). If in the previous star, HD 12932, the RV amplitude 
remained nearly constant above a certain atmospheric height, in 
y Equ we observed a continuous decrease of RV amplitude out- 
ward. The bisector amplitude-phase diagrams are fairly similar 
in both stars, showing the same constancy in the phase measured 
across Nd spectral lines and a phase shift between lines (Fig.l20l- 
Online only). This phase shift is larger in y Equ. Again, Y n lines 
have very small RV amplitudes, but their pulsation phases are 
the same as in the strongest Ndm and Prm lines. Extra broad- 
ening of the REE lines in y Equ was discussed by Kochukhov 
& Ryabchikova (2001a), who showed that y macro = 10 kms~' is 
needed for fitting strong Ndm and Prm lines. 

6.8. HD 19918 

Pulsational characteristics of HD 19918 (Fig. [22]) are similar to 
both HD 12932 and, in particular, to y Equ. An even more rapid 
decrease in the RV amplitude is observed in the atmosphere 
above a certain layer. The three stars show the largest RV ampli- 
tudes and the largest variation in the bisector velocity amplitudes 
from the line centre to the line wings (Fig. [2]]- Online only). In 



6.9. HD 9289 

The amplitude-phase diagrams for the centre-of-gravity mea- 
surements in the spectrum of HD 9289 are shown in Fig. [8] 
Because of the low S/N and significant rotation compared to 
other stars, we did not perform bisector measurements. The 
amplitude-phase diagrams for this star are similar to those for 
HD 12932. 




Fig. 8. The same as in Fig. 0but for HD 9289. 



6.10. HD 137949 (33 Lib) 

33 Lib shows the most complex pulsational behaviour among 
all roAp stars. Mkrtichian et al. (2003) find nearly anti-phase 
pulsations of Nd n and Nd m lines, which they attribute to the 
presence of pulsation node high in the atmosphere of 33 Lib. 
This was confirmed by Kurtz at al. (2005b ), who also find that in 
some REE lines the main frequency, corresponding to 8.27 min, 
and its harmonic have almost equal RV amplitudes. The high ac- 
curacy of the present observations allows us to study pulsational 
characteristics of the 33 Lib atmosphere in more detail. Figure 
[9] shows RV amplitudes (top) and phases (bottom) measured for 
the main period as a function of the central intensity of spec- 
tral lines. The NLTE Nd stratification analysis by Mashonkina 
et al. (2005 ) for HD 24712 shows that central intensities of both 
Ndn and Ndm lines are roughly proportional to there centre- 
of-gravity depth formation. Taking into account that the REE 
anomaly on which Nd stratification was based is the common 
feature for all roAp stars (Ryabchikova et al. 120011 12.0041 ). we 
can expect similar REE gradients in the atmosphere of 33 Lib, 
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Fig. 9. Amplitude (top) and phase (bottom) dependence on the central line intensity in the atmosphere of 33 Lib for the main 
frequency. Left panels represent measurements in Nd n- m lines, while right panels show the same for other elements with lines in 
the second and third ionisation states. 



hence, the similar dependence of the central intensities on the 
optical depth. 

The phase jump is defined very well by Nd lines and is also 
indicated by the measurements of Pr, Tb, Dy lines. Our results 
show that the phase jump (radial node) does not separate the 
formation regions of the lines in consecutive ionisation states as 
is claimed by Kurtz at al. ( 2005b ). Nd m, Pr m, and Th m lines are 
observed from both sides of the jump. Moreover, a position of 
the phase jump relative to line depth provides some evidence that 
abundance distributions for at least Pr, Nd, Tb, Dy are similar, 
hence one may expect abundance jumps at approximately the 
same optical depths in 33 Lib atmosphere. 

Bisector measurements support the results obtained for the 
main frequency from centre-of-gravity measurements and allow 
us to carry out a more detailed study of the pulsation properties 
in 33 Lib. Figure [10] shows bisector amplitudes (top) and phases 
(bottom) measured in a few representative spectral lines for the 
main period 8.27 min (left panels) and its harmonic (right pan- 
els). 

We have to emphasise here a major difficulty in the bisec- 
tor measurements for the lines splitted in strong magnetic field. 
Figure QT| shows intrinsic theoretical profiles of the lines from 



Fig. [9] calculated with the abundances that reproduce the ob- 
served line intensity. Only the Nd n 5320 A line does not exhibit 
resolved Zeeman components in the intrinsic profile, therefore 
this line may provide the 'purest' pulsation distribution across 
the line profile. Unfortunately, no strong Ndm or Prin line has 
negligible Zeeman splitting. The best and strongest Prin 5300 A 
and Nd in 5294 A lines both have significant splitting and rather 
different Zeeman patterns. Moreover, the Ndm 5294 A line is 
blended in both wings, and these blends may distort the observed 
pulsation effects. In particular, it concerns the RV amplitudes 
and, to a lesser extent, pulsation phases. Therefore we chose an- 
other Nd in 6327 A line as a representative of Nd in strong lines. 
Note that spectral lines with similar Zeeman patterns and simi- 
lar intensities (Nd m 5677, 5802, 585 1 A) have the same bisector 
velocity amplitude and phase distributions. Pulsation phases are 
still similar for the lines with similar Zeeman patterns but with 
different intensities (Nd m 6327 A and Prm 5300 A). Therefore, 
here we discuss mainly pulsation phase distribution in the atmo- 
sphere of 33 Lib. 

Across the profiles of all strong Prm and Ndm lines and the 
Ndn 5320 A line, one phase jump is observed for the main pe- 




Fig. 10. Amplitude-phase diagrams for the pulsation bisector measurements in 33 Lib. Amplitudes and phases are shown for the 
main period 8.27 min (left panels) and its first harmonic (right panels). 



riod and two phase jumps are detected for the first harmonic. 
Kurtz et al. (1200 5b) find only one phase jump for the first har- 
monic using bisector measurements of the Ndm6145 A line. 
The best pulsation distribution is defined by the bisector mea- 
surements in the Pnn5300 A line, where it is clearly seen that 
the position of the phase jumps corresponds to nearly zero RV 
amplitudes, indicating two radial nodes. 

Due to the diversity of the pulsational characteristics, it is 
difficult to compare pulsation phases for the lines of differ- 
ent species. However, a phase distribution for the main period 
(Fig-13 bottom left panel) suggests the existence of the 7r-radian 
jump between Ndn and Ndin lines, which was found earlier by 
Mkrtichian et al. (2003), and another 7r-radian jump between the 
Eu ii 6645 A and Dy in 5730 A lines. Comprehensive analysis of 
the variability of many REE lines in 33 Lib allowed us to ob- 
tain a refined picture of the pulsation node in high atmospheric 
layers. In particular, we show (see Fig. [9] upper panel) that Ndn 
and Nd in lines do not simply show a 7r-radian phase difference, 
but a continuous phase trend exists in the RV curves of Nd m 
lines, with some of the doubly ionised Nd lines (e.g., 5286 and 
6690 A) pulsating with phases typical of Nd n. This appears to 
be the first demonstration of the existence of a pulsation node 
within the formation region of the lines belonging to the same 
REE ion. Another remarkable observation that we made is that 
harmonic RV variations are the strongest in the lines that form 



close to the position of the node. Thus, the mechanism exciting 
the first harmonic is probably directly related to the presence of 
a radial node in the upper atmosphere of 33 Lib. 

The phases of the Ha core lie between the phases of the Eu n 
line and the wing of the Dy in 5730 A line. Interestingly, the very 
core of Ha has the highest RV amplitude, 500 m s _1 , and the last 
measured point in the Dym line wing has very high RV am- 
plitude, too, although this measurement is not accurate enough. 
Nevertheless, this gives us a hint of the line formation depths 
in the atmosphere of 33 Lib. Without detailed analysis of the 
chemical structure of the atmosphere, it is difficult to suggest a 
pulsational scenario for the star, but it is evident that we observe 
several waves over the whole atmosphere. 

The phase distribution across the Prm and Ndm spectral 
lines allows us to estimate the pulsation wave speed. Assuming 
that we have one full harmonic wave over the Nd in and Pr m 
lines and considering the geometric size of the line-formation re- 
gion in the stratified case (Mashonkina et al. 2005 1 Ryabchikova 
et al. l2006bi l. we estimated the pulsation wave speed as less than 
6 kms 1 for the harmonic period 4.136 min. It is lower than the 
sound speed in the 33 Lib atmosphere in adiabatic approxima- 
tion and is similar to the pulsation wave speed in another roAp 
star, HD 24712 (Ryabchikova et al. [20071 . The magnetic field, 
which is stronger in 33 Lib than in any other programme star, 
may be a reason for the difference in the observed distribution 



of the pulsational characteristics over the atmospheres of 33 Lib 
and HD 24712. 

What conclusions can be made from these results? 
Calculations of linear and non-linear radial pulsations (for ex- 
ample, Fadeev & Fokin 1985 ) in adiabatic approximation show 
that the wave speed may differ from the sound speed because of 
the reflection conditions and of the finite amplitude. Therefore, 
our results cannot give any proof of the non-adiabaticity of the 
pulsations in roAp stars. However, that the pulsation wave speed 
is close to the sound speed strongly supports an acoustic nature 
for the pulsations and an absence of the shock waves. 

While no extra broadening is required to fit the Fe-peak lines, 
3-4 kms~' is necessary for Yn, Eun, Gdn, and Ern lines. For 
Ndn and weak Ndm and for Prin lines, one needs 6-8 kms -1 , 
and, finally, for the strong Nd in and Pr m lines we have to intro- 
duce V maC ro = 10-12 km s _1 to fit the observed profiles. 

Rather high amplitude pulsations were detected in both Th in 
lines with the phases corresponding to those for Tb m lines. The 
RV amplitudes of the main frequency derived for the lines of 
Na, Mg, Si, Ca, Sc Cr, Ti, Fe, Ni lie in the 30-50 ms" 1 range, 
while the amplitude of the harmonic is below the detection limit 
of ~10 ms -1 . The Ba lines also reveal pulsations at the main 
frequency, but with higher RV amplitudes (70-100 ms~'). The 
Y ii lines have similar amplitudes, which seem to grow with the 
line intensity. The Y n lines also show detectable harmonic am- 
plitudes. All these lines show pulsation phases in the range of 
pulsation phases for the singly ionised REEs. 
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Fig. 11. Zeeman structure for a set of REE spectral lines in 
33 Lib. 



7. Summary and discussion 

We have analysed spectroscopic pulsational variability of ten 
roAp stars in detail. For each object, several hundred spectral 
lines were measured and the resulting time series of radial ve- 
locity and bisector variation of selected lines were interpreted 
with the help of the least-square fitting technique. We confirm 
results of the previous studies, which suggest that non-radial os- 
cillations are primarily detectable in the lines of heavy elements, 
especially in the doubly ionised REE lines. This phenomenon is 
attributed to the extreme chemical stratification present in the at- 
mospheres of cool Ap stars. Due to the segregation of elements 
under the influence of radiative diffusion and due to lack of mix- 
ing, a thin layer enriched in heavy elements is created in the up- 



per atmosphere, where rather weak roAp non-radial oscillations 
attain significant amplitudes. 

Stellar atmosphere modelling that includes stratification, 
NLTE, and magnetic field effects is an extremely complex prob- 
lem. Up to now only the first preliminary models of different 
element distributions over stellar atmosphere based on the spec- 
troscopic observations were calculated for a few cool Ap stars. 
Self-consistent diffusion calculations (LeBlanc & Monin 2004) 
provide the next step in our understanding of the Ap star phe- 
nomenon. However, with the diffusion calculations for 39 ele- 
ments from He to La, LeBlanc & Monin's models do not include 
REE elements, which show the most outstanding stratification 
signatures in roAp atmospheres. The predicted stratification of 
La is far from the observed one because a small number of the 
ionisation states (two) is taken into account and NLTE effects are 
not considered. Empirically derived Nd abundance distribution 
shows directly that NLTE effects play a crucial role in the REE 
stratification studies (Mashonkina et al. 2005). Based on the em- 
pirical element distribution in the atmospheres of two roAp stars, 
HD 24712 and y Equ, we conclude that pulsation amplitudes 
and, in particular, phases derived from the lines of different ele- 
ments correlate with the optical depth. Therefore, in the absence 
of the proper atmospheric models, the amplitude-phase diagrams 
proposed in the present paper are a powerful tool in the study of 
vertical structure of p-modes in the atmospheres of roAp stars. 

Our pulsation analysis of the radial velocity variations 
demonstrates similarity in the atmospheric pulsation characteris- 
tics. With the help of amplitude-phase diagram analysis we find 
that pulsation waves exhibit either a constant phase and ampli- 
tude changing with height or depth-dependence of both parame- 
ters. We interpret the former as a signature of standing pulsation 
wave and the latter as evidence of a travelling (running) wave 
in stellar atmosphere. Thus, in general, atmospheric pulsational 
fluctuations in roAp stars can be represented by a superposition 
of standing and running waves. According to our results, pul- 
sation waves in three roAp stars, HD 24712, HD 134214, and 
a Cir, with the pulsation frequency close to or below the acous- 
tic cut-off limit, have the running-wave character from the low 
atmospheric heights. In the longer period stars, standing waves 
are observed up to some of the atmospheric layers, defined by 
the formation depth of the lines of specific elements, while run- 
ning waves dominate higher in the atmosphere. 

We also find that the change from the standing to running 
character of pulsation waves depends on the effective tempera- 
ture of roAp star. Hotter stars seem to develop the running wave 
deeper in the atmosphere than the cooler stars with the same pul- 
sation periods. 

In all stars but 33 Lib, independent of the atmospheric and 
pulsation parameters, pulsation measurements reveal waves trav- 
elling through the layers defined by the same sequence of chem- 
ical species. The lowest amplitudes are observed for Eun lines, 
then pulsations propagate in the layers where Ha core, Nd, and 
Pr lines originate. Pulsation amplitude reaches maximum around 
these atmospheric heights and then decreases outwards in most 
stars. The RV extremum of the second REE ions is always ob- 
served later in time relative to the variation of singly ionised 
REEs. The largest phase shifts and amplitudes are often detected 
in Tbm and Thin lines. Pulsational variability in the latter lines 
is detected here for the first time. Similarity of the pulsation 
wave propagation signatures in the studied roAp stars suggests 
that layers enriched in different REE species are arranged in ap- 
proximately the same vertical order in all stars. For two stars, 
33 Lib and HD 19918, we find weak but definite RV variability 



in the lines of iron-peak elements, confirming previous results by 
Mkrtichian et al. (120031 1 and Kurtz et al. (2005b), respectively. 

In all stars, spectral lines with the highest RV amplitudes 
have an additional broadening varying from 4 to 10-12 kms -1 
in terms of the macroturbulent velocity required to reproduce the 
observed line profile shapes. The running wave characteristics 
usually appear in the lines with V maC ro ^10 km s . According to 
Kochukhov et al. ((2007 ), these are the lines where the usual sym- 
metric pulsational line profile variability is transformed into an 
asymmetric blue-to-red moving pattern. These authors propose 
pulsational variability of the line widths, arising from the peri- 
odic expansion and compression of the turbulent layers in the 
upper atmospheres of roAp stars as an explanation for the asym- 
metric line profile variability pattern. The present study shows 
that the position of the turbulent layer in the roAp-star atmo- 
spheres is defined by the formation depth of Pr in lines in a cooler 
part of the programme stars and by the formation depth of Nd m 
lines in hotter stars. This turbulent layer, which is probably re- 
lated to the REE abundance gradients in the upper atmosphere 
(Mashonkina et al. 2005 ), seems to be the key element in the 
depth-dependence of the spectroscopic pulsation characteristics 
in roAp stars. Abundance gradients may cause a non-standard 
temperature structure in Ap-star atmospheres. This is supported 
by an atmospheric anomaly that required to explain the pecu- 
liar core-wing transition of the Balmer lines (Kochukhov et al. 
2002). Therefore, the comprehensive investigation of the roAp- 
star atmosphere, including chemical diffusion and magnetohy- 
drodynamic modelling of the interaction of convection, pulsa- 
tion and magnetic field, is necessary to fully understand the full 
variety of the phenomena associated with a turbulent REE-rich 
cloud in cool Ap stars. 

One of the sample stars, 33 Lib, differs from all other roAp 
stars in that, at a given atmospheric heights, it shows a compara- 
ble RV amplitude of the main frequency and its first harmonic. 
Comprehensive analysis of the REE lines of several elements in 
both ionisations stages shows that harmonic oscillation emerges 
close to the position of the pulsation node located within the 
REE-rich high atmospheric layer. This follows from the obser- 
vation that the lines showing large double-wave variation are all 
located close to the minimum amplitude and to a 7r-radian phase 
jump of the main frequency. Thus, the physical mechanism giv- 
ing rise to harmonic spectroscopic variability in 33 Lib must be 
closely related to the existence of the radial node. 

Based on the observation of the two radial nodes for the har- 
monic variation across the profiles of a few strongest Nd in and 
Prm lines, we can estimate a radial wavelength and pulsation 
wave speed. The latter is below the adiabatic sound speed in 
the atmosphere of 33 Lib. The same result was obtained ear- 
lier for another roAp star, HD 24712, using a different approach 
(Ryabchikova et al. 2007 ). These results support an acoustic na- 
ture for the pulsations and reject the idea of the shock wave pro- 
posed by Shibahashi et al. ([2004;) for interpreting of the blue- 
to-red pulsation pattern, because this model requires supersonic 
pulsation motions. 

33 Lib is the only star showing definite and direct evidence of 
radial nodes in the atmosphere. But the absence of a measurable 
pulsation signal in the lines of the elements with Z < 38 in most 
of our programme stars poses a question about the existence of 
another node in the lower photospheric layers. Abundance strat- 
ification analyses of cool Ap stars (Ryabchikova et al. 2002 
2005b; Kochukhov et al. 2006 ), as well as the model atmosphere 
with self-consistent diffusion calculations (LeBlanc & Monin 
2004), show that light and iron-peak elements have a tendency 
to concentrate in the deeper atmosphere, below log T5000 = — 1, 



and to be strongly depleted in the outer atmospheric layers. 
Therefore, practically all Fe-peak lines observed in the optical 
spectral region are formed in the layers -1.5 < logrsooo < 0.0. 
The non-adiabatic theoretical pulsation model for HD 24712 
supports the existence of the nodal region near the photosphere 
(logTsooo ~ 0) (see Fig. 3 in Sachkov et al. 2006). The position 
of the radial node depends on the frequency (or the radial order) 
of the mode. As the frequency increases, the node shifts out- 
wards (Saio, private communication). According to this model, 
in the roAp stars HD 24712, HD 134214, and a Cir, the node 
is located close to, or slightly above, the continuum-forming 
layers. On the other hand, in the longer-period stars the node 
should be located below the photosphere (see Figs. 3, 4 in Saio 
& Gautschy 2004; Fig. 8 in Saio 2005). For some intermediate 
pulsation frequencies, we should detect pulsation signal in the 
region of the iron-peak line formation, but this is not observed. 
Perhaps, the pulsation amplitude in the lower atmosphere is di- 
minished to such an extent that pulsations are undetectable even 
in anti-node regions. 

A discovery of the low-amplitude pulsations in Y 11 lines in 
phase with the highest-amplitude Prm lines delineates another 
problem. In the three stars with the shortest pulsation periods, 
HD 24712, HD 134214, and a Cir, the Yn lines produce a sec- 
ondary minimum in the amplitude-phase diagrams. These lines 
have the lowest detected RV amplitudes and are out of phase by 
7r-radian with some other low-amplitude lines, for instance, Eu n. 
Does this indicate a radial node near the Eun line formation 
heights or do we see a double-wave structure in the amplitude- 
phase diagram depending on the chemical structure of stellar at- 
mospheres? 

The amplitude-phase diagrams derived for as many elements 
and spectral lines as possible are proven to be an extremely pow- 
erful tool for investigating the pulsation properties of roAp atmo- 
spheres. Detailed inferences about the vertical mode structure 
obtained in our study call for in-depth theoretical investigation 
of the propagation of pulsation waves in magnetically dominated 
and chemically stratified atmospheres. The most important and 
difficult theoretical challenge is to recognise and model the phys- 
ical processes that are responsible for remarkable observation of 
pulsation wave transformation in the REE-rich layer. Why does 
the dominant character of pulsational perturbation changes from 
standing to running wave? What causes pulsation amplitude to 
diminish about certain atmospheric height? Why does isotropic 
turbulence increase dramatically in this layer? And finally, what 
is the origin of the complex bisector behaviour observed in sev- 
eral roAp stars? Substantial theoretical developments are needed 
to resolve these issues. 
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Fig. 15. The same as in Fig.|5]but for HD 24712. Fig.21. The same as in Fig.|5]but for HD 19918. 
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Fig. 18. The same as in Fig.[3]but for HD 12932. 
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Fig. 20. The same as in Fig.|5]but for HD 201601. 
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Fig. 12. Bisector measurements as a function of the line residual intensity in HD 101065. Left panel represents RV amplitudes, ri| 
panel shows phases calculated with the main period from Table 01. 
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Fig. 13. The same as in Fig.|4]but for HD 122970. 




Fig. 14. The same as in Fig.g]but for HD 24712. 
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Fig. 17. The same as in Fig.fflbut for HD 12932. 
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Fig. 19. The same as in Fig.|4]but for y Equ. 




Fig. 22. The same as in Fig.g]but for HD 19918. 



Table 3. Pulsation amplitudes (first line) and phases (second line) calculated with the period highlighted with italics in Table[2] Errors in the last 
digits are given in paranthesis. When no pulsation signal was detected, we give the formal amplitude solution for the fixed period without phase 
information 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


Fe I - Fe II 


5405.775 


7(5) 


7(4) 


14(4) 


11(4) 


25 (18) 


37(15) 
0.607 (67) 


17 (19) 


9(3) 


11(20) 


17( 2) 
0.686(21) 


5410.910 


12(4) 




11(4) 


6(3) 


16(14) 


52(18) 
0.650 (56) 


60 (24) 


6(3) 


23(12) 


14(1) 
0.632(22) 


5415.199 


8(4) 


21(7) 


2(4) 


9(3) 


18(13) 


21(9) 
0.624(70) 


14 (14) 


5(4) 


48(23) 


13(1) 
0.396(24) 


5424.068 


4(4) 




6(4) 


2(3) 


13( 9) 


29( 8) 
0.640(46) 


31 (18) 


2(5) 


38(24) 


25(2) 
0.685(14) 


5434.524 


20(7) 


51 (21) 


9(7) 


4(3) 


7(15) 


47 (19) 
0.665 (65) 


38 (25) 


11(6) 


95 (44) 


27(2) 
0.669 (15) 


5576.089 


15 (10) 




27 (10) 


13(4) 


35(15) 


100 (43) 
0.778 (69) 


33 (27) 


11(8) 


126(41) 


22(2) 
0.642 (16) 


5414.073 


18(13) 


11(14) 


5( 8) 


17(8) 


30(17) 


34(9) 
0.707(41) 


50 (31) 


6(6) 


24(23) 


17( 2) 
0.552(26) 


5425.257 


36(12) 




24(8) 


11(6) 


41(12) 


35(17) 
0.545 (77) 


65 (23) 


23 (11) 


48(34) 


25(3) 
0.593 (21) 


6247.557 




13(11) 


1(5) 


34(14) 


63(13) 


29( 7) 
0.722(43) 




44( 9) 


42(43) 


35(3) 
0.646(15) 


YII 


3549.005 


74(10) 
0.722 (23) 




















4883.684 














102 (30) 
0.657 (48) 








4982.129 


19(8) 


69(18) 
0.824 (42) 








86 (55) 










5087.416 


53(5) 




120 ( 4) 


73(5) 


132(19) 


89 (23) 


114(36) 


109 ( 8) 


133 (30) 


61 (2) 




0.633 (17) 




0.893 ( 6) 


1.141 (11) 


1.182(24) 


1.102(42) 


0.700 (51) 


0.596 (12) 


0.623 (36) 


0.609 ( 8) 


5119.112 


17(7) 




27(18) 


38 (10) 
1.031 (43) 














5123.211 






82(6) 
0.860 (13) 
















5200.406 






86(8) 
0.917 (16) 






62 (18) 
1.171 (49) 










5205.724 






107 ( 5) 
0.890 ( 8) 






57(15) 
0.893 (44) 










5402.774 


29(6) 




92(13) 
0.904 (23) 


53(5) 
1.025 (18) 






211 (59) 
0.568 (45) 






78(4) 
0.688 ( 9) 


5473.388 






68 (16) 
0.831 (38) 
















5509.895 


31(7) 


36(6) 


83(6) 


19(5) 










66 (27) 


49(3) 




0.556 (37) 


1.012(27) 


0.882 (12) 


1.021 (42) 










0.489 (66) 


0.610(10) 


5544.611 








34(9) 
1.285 (42) 






113 (57) 
0.833 (80) 








5546.009 












118(24) 
0.737 (33) 








94(4) 
0.723 ( 7) 


5662.925 


37(6) 




131(6) 


71(6) 




48 (16) 


107 (25) 


121 ( 3) 


49(16) 


98(3) 




0.568 (26) 




0.901 ( 8) 


1.101 (14) 




0.995 (54) 


0.615 (38) 


0.604 ( 4) 


0.584 (52) 


0.677 ( 6) 


5728.890 






140 (30) 
0.896 (34) 


42 (14) 
1.054 (55) 








56 (21) 




104 ( 7) 
0.640(11) 


6613.733 






98 (20) 
0.903 (33) 










77 (12) 
0.589 (25) 






6795.414 






158(19) 
0.858 (19) 


81 (14) 
0.880 (28) 








59 (14) 
0.477 (39) 




73(8) 
0.657 ( 1) 


La II 


4804.039 


110(9) 
0.296(13) 












373 ( 94) 
0.529 (41) 








4921.776 














500 ( 53) 
0.530 (17) 








4986.819 


78 (9) 
0.325 (19) 




















5114.559 






257 (15) 
0.837 ( 9) 


265 ( 6) 
0.765 ( 4) 














5122.988 






217(12) 
0.854 ( 9) 
















5157.433 






183 (23) 
0.836 (20) 


235 ( 9) 
0.734 ( 7) 














5163.611 






226 (27) 
0.839(19) 


207 (11) 
0.697 ( 9) 














5259.379 






227 (20) 
0.841 (14) 


196 ( 9) 
0.686 ( 8) 














5290.818 








188 ( 9) 
0.733 ( 8) 






448 ( 73) 
0.572 (26) 








5303.528 


100 ( 8) 
0.274(13) 




194(14) 
0.849 (12) 


219 ( 9) 
0.744 ( 7) 














5340.655 








188(11) 
0.739 (10) 














5377.052 








260 ( 8) 
0.744 ( 5) 















Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


5458.697 








315 (19) 
0.796 ( 9) 














5482.268 








221 (18) 
0.727 (13) 














5566.925 




121(6) 
0.849 ( 9) 


















5797.565 






202(16) 
0.820(13) 
















5805.773 




126(5) 
0.833 ( 7) 


239 (20) 
0.838 (14) 


224(12) 
0.756 ( 9) 






396(134) 
0.463 (55) 






111(9) 
0.662(13) 


5808.313 




66(6) 
0.834(14) 


215 (22) 
0.827 (17) 


261 (15) 
0.731 (9) 












174 (10) 
0.630 ( 9) 


5863.691 








235 (13) 
0.740 ( 9) 














5880.633 




94(6) 
0.814(10) 




213(15) 
0.733 (11) 














5936.210 








273 (15) 
0.718(9) 












154(9) 
0.704(10) 


6100.374 






166 (30) 
0.842 (29) 


183(13) 
0.715 (11) 














6129.556 




56(6) 
0.815(17) 


















6146.523 








180 (27) 
0.745 (24) 














6262.287 


81 (12) 
0.257 (23) 


















198 (12) 
0.672(10) 


6296.079 


105 (11) 
0.317(17) 




















6310.912 




46(5) 
0.789 (20) 


















6320.376 


78(11) 
0.288 (23) 




195(18) 
0.858 (15) 


213(8) 
0.730 ( 6) 












114(7) 
0.634 (10) 


6358.105 




88(4) 
0.809 ( 8) 


















6390.477 








205(11) 
0.754 ( 9) 














6399.030 


73(18) 
0.223 (41) 


74(9) 
0.821 (21) 




295 (25) 
0.739 (14) 












177 (16) 
0.643 (15) 


6642.755 








274(17) 
0.723 (10) 












111 (7) 
0.660 (10) 


6808.885 








251 (56) 
0.702 (35) 














PrII 


5070.023 




175 (11) 
0.863 (10) 
















88(7) 
0.892(13) 


5110.760 


200(12) 
0.446(10) 




















5129.533 


243(11) 
0.477 ( 8) 




















5135.140 


313(16) 
0.454 ( 9) 




745 (60) 
0.865 (13) 


332 (37) 
0.980(18) 




329 (136) 
0.940 (65) 


734(170) 
0.570 (37) 








5207.895 




310(7) 
0.867 ( 4) 


















5219.045 






555 (40) 
0.892(11) 
















5292.619 


195 (18) 
0.400(14) 


















131(6) 
0.807 ( 8) 


5321.817 




149 ( 5) 
0.852 ( 6) 


















5322.772 


368 (18) 


493 ( 6) 


872 (35) 


480 (23) 


392 (101) 


390 (167) 


704(191) 










0.468 ( 8) 


0.895 ( 2) 


0.866 ( 6) 


0.969 ( 8) 


1.028 (39) 


0.859 (67) 


0.605 (44) 








5432.038 




2215 (10) 
0.818(8) 


















5519.388 




308 (10) 
0.836 ( 6) 


















5605.642 




337 ( 6) 
0.843 ( 3) 




276 (53) 
0.881 (31) 














5608.931 




359 ( 5) 
0.821 (7) 
















69 (12) 
0.771 (28) 


5681.901 


159 (44) 
0.355 (44) 






346 (72) 
0.807 (33) 












325 (29) 
0.811 (14) 


5786.166 




460 (17) 
0.871 (6) 


















5879.260 




382 ( 9) 
0.856 ( 4) 


















5939.899 




436 (10) 
0.872 ( 4) 


















6017.767 


170 (28) 
0.385 (26) 




















6148.243 




378 (18) 
0.841 (8) 


















6165.891 


251 (14) 
0.456 ( 9) 


404 ( 8) 
0.860 ( 3) 


711 (40) 
0.865 ( 9) 


267 (38) 
1.003 (23) 












186 ( 8) 
0.812(7) 



Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


6197.452 




392 (26) 






















0.829 (11) 


















6255.098 




304 (15) 






















0.816(8) 


















6656.834 


254 (26) 






















0.419(17) 




















Prill 


4910.823 


248 (11) 






















0.543 ( 7) 




















5284.693 


259 (1 1) 


427 ( 4) 


600 (14) 


380 (12 


1043 ( 33) 


925 ( 40) 


601 ( 27) 


496 ( 8) 


504 ( 21) 


71 (4) 




0.642 ( 7) 


0.934 ( 2) 


0.887 ( 4) 


1.221 ( 5) 


1.211 ( 5) 


1.042 ( 7) 


0.804 ( 7) 


0.692 ( 3) 


0.633 ( 7) 


0.263 (10) 


5299.993 


246 (10) 


460 ( 5) 


596 (14) 


365 (13) 


1192 ( 32) 


1 185 ( 52) 


95 1 ( 47) 


491 ( 15) 


724 ( 30) 






0.641 ( 7) 


0.946 ( 2) 


0.888 ( 4) 


1.209 ( 6) 


1.189 ( 4) 


1.055 ( 7) 


0.814 ( 8) 


0.684 ( 5) 


0.677 ( 7) 




5550.593 












325 (102) 






















0.958 ( 50) 










5765.260 




589 (15) 


1005 (60) 


300 (40) 


847 ( 83) 


943 (107) 


578 (133) 


420 ( 66) 


906 (173) 


85 (15) 






0.895 ( 4) 


0.871 ( 9) 


1.072 (22) 


1.130 ( 15) 


1.003 ( 18) 


0.640 ( 36) 


0.587 ( 25) 


0.581 ( 31) 


0.894 (28) 


5844.408 


299 (15) 


398 ( 7) 


753 (30) 


379 (19) 


746 ( 79) 


911 ( 57) 


599 ( 94) 


414 ( 37) 


490 ( 54) 


11(5) 




0.576 ( 8) 


0.918 ( 3) 


0.881 ( 6) 


1.120 ( 8) 


1.173 ( 17) 


1.076 ( 10) 


0.742 ( 24) 


0.677 ( 15) 


0.591 ( 18) 


0.473 (76) 


5956.050 


211 ( 8) 












589 ( 41) 










0.616 ( 7) 












0.758 ( 11) 








5998.968 


329 (15) 












614 (122) 










0.587 ( 8) 












0.742 ( 31) 








6053.004 


307 (12) 






















0.583 ( 7) 




















6090.010 


265 (12) 




685 (18) 




1072 ( 31) 










46(2) 




0.631 ( 8) 




0.885 ( 4) 




1.168 ( 5) 










0.316 ( 9) 


6160.233 


251 (10) 




655 (19) 


376 (17) 


1408 ( 39) 


1241 ( 34) 




441 ( 14) 


752 ( 43) 


53 ( 5) 




0.624 ( 6) 




0.888 ( 5) 


1.169 ( 7) 


1.161 ( 4) 


1.025 ( 4) 




0.666 ( 5) 


0.679 ( 9) 


0.225 (16) 


6195.621 


272(12) 






335 (11) 


1185 ( 32) 


1036( 37) 




481 ( 13) 


734(91) 


61 ( 3) 




ft f\0f\ ( 1\ 






1 1 77 ( 5*1 


1 155 ( 4) 


1 016 1 f,\ 
i.uio t O) 




676 C A\ 


fi 79^ ( 9fYt 
U. 1 Ld \ ZAj) 


\).z.o I \ y ) 


6429.265 




352 ( 9) 






















0.902 ( 4) 


















6500.044 


244(12) 




638 (37) 




720 ( 67) 


944 ( 53) 












0.527 ( 8) 




0.888 (9) 




1.098 ( 15) 


0.993 ( 9) 










6578.901 




531(8) 


1071 (46) 




889 (123) 


1175 (141) 






792(139) 








0.893 ( 3) 


0.886(7) 




1.092(22) 


1.025 ( 19) 






0.561 (28) 




6692.247 


317 (27) 


418 (15) 


862 (89) 


244 (55) 


495 (137) 


687 (134) 




689 (234) 




94(11) 




0.501 (14) 


0.884 ( 6) 


0.855 (17) 


1.085 (36) 


1.080(44) 


0.925 ( 31) 




0.655 ( 53) 




0.870(19) 


6706.704 


359(15) 


530 (10) 


817 (33 


372 (23) 


1027 ( 93) 


1353 ( 79) 




553 ( 52) 


773 (118) 


49(8) 




0.569 ( 7) 


0.918(3) 


0.882(6) 


1.135(10) 


1.132 ( 14) 


1.017 ( 9) 




0.630 ( 15) 


0.608 ( 25) 


0.929 (28) 


6866.794 




527 ( 4) 


806 (27) 


405 (13) 








816 ( 18) 




66(5) 






0.927 ( 1) 


0.885 (5) 


1.158(5) 








0.673 ( 4) 




0.184(13) 


Nd II 


3375.244 


182 (23) 






















0.365 (20) 




















4706.543 














858 ( 56) 






















0.595 ( 11) 








4797.153 


205 ( 9) 






















0.368 ( 7) 




















4799.420 


214 ( 9) 






















0.346 ( 7) 




















4811.342 


274 ( 9) 












1087 ( 59) 










0.393 ( 5) 












0.567 ( 9) 








4947.020 














1302(117) 






















0.548 ( 14) 








4959.119 


254 ( 9) 






















0.383 ( 6) 




















4987.161 






422 (28) 






















0.840(11) 
















4989.950 


280 ( 8) 




519 (32) 






1157 ( 92) 


1260(92) 










0.363 ( 5) 




0.880 (10) 






0.838 ( 13) 


0.582 ( 12) 








5027.847 




239 ( 7) 






















0.841 (5) 


















5033.507 


257 (11) 


221 (5) 


548 (40) 


528 (19) 


825 (82) 




1037 (103) 


659 (63) 


399 (81) 


203 (10) 




0.351 (7) 


0.851 (4) 


0.822(12) 


0.866 ( 6) 


1.000(16) 




0.563 ( 16) 


0.362(15) 


0.604 (33) 


0.733 ( 8) 


5077.154 


199 ( 9) 


167 ( 4) 


571 (42) 


475 (25) 






1226 (125) 






137 ( 8) 




0.326 ( 7) 


0.841 (4) 


0.872(12) 


0.842 ( 8) 






0.565 ( 17) 






0.735 (10) 


5089.832 


214(8) 




437 (18) 








834(105) 






72(3) 




0.358 ( 6) 




0.843 ( 7) 








0.563 (21) 






0.747 ( 8) 


5092.794 


255 ( 9) 




417(17) 


495 (13) 


992 (52) 


895 ( 58) 


943 ( 80) 






59(4) 




0.378 ( 6) 




0.840 ( 6) 


0.885 (4) 


0.996 ( 8) 


0.822 ( 10) 


0.566 ( 14) 






0.705 (11) 


5096.516 


236 ( 9) 


239 ( 2) 


593 (28) 


453 (19) 


500 (68) 


440 ( 80) 








117(6) 




0.328 ( 6) 


0.842 ( 2) 


0.854 ( 8) 


0.857 (7) 


0.944(21) 


0.817(28) 








0.770 ( 9) 


5105.226 






439 (22) 














193 ( 7) 








0.838 ( 8) 














0.768 ( 6) 


5130.586 


261 ( 8) 


208 ( 2) 


448 (13) 


499 (10) 


976 (40) 


1131 ( 56) 


1145 ( 60) 


773 (17) 


415 (36) 






0.398 ( 5) 


0.855 ( 2) 


0.858 (5) 


0.915(3) 


1.015(6) 


0.843 ( 8) 


0.593 ( 8) 


0.362 ( 4) 


0.576(14) 




5132.328 


279(11) 




477 (23) 


522(15) 


707 (48) 


898 ( 80) 








117(3) 




0.372 ( 6) 




0.862 ( 8) 


0.887 (5) 


0.993 (10) 


0.839 ( 14) 








0.778 ( 5) 


5143.337 


237 (11) 


234 ( 3) 


599 (33) 


404 (23) 


502 (72) 




1012(171) 









Table 3. Continued. 



Wavelength 
(A) 


Amplitude (m s 1 ) / Phase (period fraction) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 




0.332 ( 7) 


0.827 ( 3) 


0.838 (9) 


0.850(9) 


0.971 (22) 




0.538 ( 27) 








5152.62 






416(55) 


277 (56) 




















0.835 (21) 


0.799 (33) 














5165.129 


214(8) 


205 ( 3) 


374(18) 


443(11) 


618(34) 


676 ( 48) 








44(4) 




0.356 ( 6) 


0.857 ( 3) 


0.857 ( 8) 


0.864 ( 4) 


0.981 (9) 


0.839 ( 11) 








0.750 (14) 


5167.919 




244 ( 5) 






















0.823 ( 3) 


















5176.781 


215(11) 




436 (28) 






663 ( 81) 












0.330 ( 8) 




0.849 (10) 






0.799 ( 18) 










5181.169 




234 ( 3) 
















71 (4) 






0.847 ( 2) 
















0.741 (10) 


5182.587 


251 (11) 










714(72) 


745 ( 54) 






158 ( 7) 




0.336 ( 7) 










0.811 ( 16) 


0.529 ( 12) 






0.758 ( 7) 


5192.614 




197 ( 3) 






















0.849 ( 3) 


















5212.361 




291 ( 4) 


440(18) 














91 (4) 






0.843 ( 3) 


0.845 ( 6) 














0.766 ( 7) 


5255.506 


263 ( 9) 


292 ( 4) 


440(13) 














59(4) 




0.391 (5) 


0.846 ( 2) 


0.844 ( 5) 














0.762(11) 


5276.869 


219(8) 




516(21) 


489 (13) 


792 (62) 




971 ( 64) 




314(53) 


120 ( 5) 




0.358 ( 6) 




0.844 ( 7) 


0.862 ( 4) 


0.959(12) 




0.573 ( 11) 




0.611 (27) 


0.787 ( 8) 


5287.133 




239 ( 4) 






















0.834 ( 3) 


















5293.163 






387 (15) 










858 (28) 












0.855 ( 6) 










0.391 (5) 






5306.460 




228 ( 3) 


544 (26) 


446(15) 


633 (63) 
















0.849 ( 2) 


0.854 ( 8) 


0.859 ( 6) 


0.963 (16) 












5311.453 


243 ( 9) 


207 ( 2) 


493 (21) 


482(14) 


798 (58) 


1160 (102) 




785 (32) 


572(121) 


97(4) 




0.372 ( 6) 


0.839 ( 2) 


0.843 ( 7) 


0.885 (5) 


0.986(12) 


0.820 ( 14) 




0.323 ( 7) 


0.599 (34) 


0.719(8) 


5314.554 


236(14) 


210(4) 


514(48) 


391 (28) 












151 (12) 




0.313(10) 


0.839 ( 3) 


0.824(15) 


0.837 (12) 












0.776(13) 


5319.815 


262 ( 9) 


209 ( 3) 


434(13) 


519(12) 


1000 (29) 


1142(43) 


1079 ( 67) 


793 (15) 


609 ( 5) 


125 ( 5) 




0.398 ( 5) 


0.847 ( 3) 


0.846 ( 5) 


0.924 ( 4) 


1.113(5) 


0.846 ( 6) 


0.592 ( 10) 


0.356 ( 3) 


0.545 (14) 


0.793 ( 7) 


5356.967 


241 (10) 


181 (4) 


485 (23) 




1035 (75) 






858 (35) 


603 (78) 






0.359 ( 7) 


0.847 ( 4) 


0.848 ( 7) 




0.975 (11) 






0.323 ( 6) 


0.605 (21) 




5361.467 


224 ( 7) 




412(35) 










416 (18) 


341 (48) 






0.383 ( 5) 




0.830(14) 










0.319(7) 


0.493 (23) 




5371.927 




242 ( 4) 


459 (23) 


469 (14) 


701 (38) 


706 ( 45) 


1188 ( 99) 












0.839 ( 3) 


0.852( 8) 


0.891 (5) 


0.997 ( 9) 


0.868 ( 10) 


0.552 ( 14) 








5385.888 




243 ( 4) 


542 (40) 


463 (21) 


654 (69) 


641 (76) 




577 (55) 










0.830( 3) 


0.845 (12) 


0.860 ( 7) 


0.934 (17) 


0.840 ( 18) 




0.309 (15) 






5399.099 




217(6) 


565 (39) 


373 (28) 






1043 (142) 












0.826 ( 4) 


0.847 (11) 


0.817(12) 






0.495(21) 








5416.374 




222 ( 4) 


560 (28) 


467 (14) 




852 ( 80) 


943 ( 89) 


592 (61) 




112(7) 






0.854 ( 3) 


0.837 ( 8) 


0.850 ( 5) 




0.835 ( 14) 


0.550 ( 15) 


0.288 (17) 




0.710(10) 


5430.780 




244 ( 6) 


441 (67) 




















0.830 ( 4) 


0.865 (24) 
















5431.516 






422 (36) 


418(15) 


615(61) 


762 ( 96) 


1285 (180) 


718 (49) 




168 (10) 








0.864(14) 


0.881 (6) 


0.989(16) 


0.824 ( 20) 


0.534 ( 23) 


0.314(11) 




0.733 (10) 


5442.264 




201 ( 3) 


447 (19) 




















0.847 ( 2) 


0.842 ( 7) 
















5456.553 














859 (100) 




313(85) 


103 ( 6) 
















0.524 ( 19) 




0.538 (44) 


0.749 (11) 


5474.731 




198 ( 4) 


535 (38) 


417 (20) 


334 (57) 
















0.847 ( 4) 


0.850(11) 


0.843 ( 8) 


0.942 (28) 












5485.696 


230 ( 8) 


226 ( 6) 


483 (27) 


484(18) 


741 (65) 


1013(81) 




705 (24) 


238 (71) 


109 ( 4) 




0.379 ( 6) 


0.839 ( 4) 


0.860 ( 9) 


0.891 (6) 


0.975 (14) 


0.837 ( 13) 




0.333 (6) 


0.539 (48) 


0.736 ( 6) 


5491.390 




162 ( 6) 






















0.806 ( 6) 


















5533.820 


245 (13) 




500 (38) 


492 (32) 


589 (99) 


1215 (170) 


1126 (121) 


921 (85) 




128(7) 




0.345 ( 9) 




0.816(12) 


0.863 (11) 


0.974 (31) 


0.781 (21) 


0.503 ( 17) 


0.287 (15) 




0.771 (9) 


5557.40 






505 (35) 














164 ( 8) 








0.829 (11) 














0.782( 8) 


5569.946 




263 ( 4) 






















0.835 ( 3) 


















5618.994 






545 (42) 


413 (20) 


346 (91) 


638 ( 96) 
















0.857 (12) 


0.846 ( 8) 


0.932 (43) 


0.846 ( 24) 










5659.775 






401 (19) 






















0.859 ( 7) 
















5687.957 




255 ( 4) 






















0.839 ( 3) 


















5702.238 




202 ( 3) 
















128 ( 6) 






0.829 ( 3) 
















0.748 ( 7) 


5707.363 




221 (7) 






















0.809 ( 5) 


















5718.180 




212(3) 


427 (23) 














69(3) 






0.815(3) 


0.854 ( 9) 














0.698 ( 8) 


5726.825 








374(14) 


527 (56) 




















0.850 ( 6) 


0.986 (17) 












5734.550 


245 (21) 


259 (10) 


386 (35) 


363 (24) 




451 (117) 




378 (92) 




152(12) 




0.296(14) 


0.833 ( 7) 


0.830(15) 


0.845 (11) 




0.798 ( 40) 




0.259 (39) 




0.750(13) 


5740.858 






472 (36) 










853 (122) 












0.844(12) 










0.320 (23) 







Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


5742.087 




239 ( 5) 
0.835 ( 4) 


503 (34) 
0.857 (11) 
















5743.193 




239 ( 6) 
0.827 ( 4) 


392 (53) 
0.814 (22) 


459 (38) 
0.821 (13) 












116(8) 
0.782(11) 


5744.777 
















508 (59) 
0.266 (19) 


477(156) 
0.432 (53) 


139 ( 9) 
0.749 (10) 


5761.694 




253 ( 7) 
0.825 ( 4) 




346 (20) 
0.849 (10) 








487 (88) 
0.275 (29) 




137 (10) 
0.748 (12) 


5764.212 




221 (6) 
0.819(4) 


328 (41) 
0.857 (19) 
















5804.004 








404(11) 
0.865 ( 5) 


573 (52) 
0.938 (14) 


608 ( 46) 
0.825 ( 12) 






321 (49) 
0.512(25) 


108 ( 4) 
0.759 ( 7) 


5811.572 




217(5) 
0.826 ( 4) 


588 (26) 
0.848 ( 7) 


456(17) 
0.860 ( 6) 


768 (82) 
0.959 (17) 












5825.857 




215(4) 
0.836 ( 4) 


479 (26) 
0.846 ( 9) 


457 (16) 
0.851 (6) 










286 (85) 
0.592 (48) 


119(5) 
0.768 ( 7) 


5842.366 








338 (16) 
0.884( 8) 












108 (11) 
0.837 (17) 


5865.027 




252 ( 7) 
0.825 ( 5) 


462 (42) 
0.859 (15) 


344 (33) 
0.878 (16) 












142(11) 
0.750(13) 


5891.506 




226 ( 3) 
0.831 (3) 


417 (28) 
0.844(10) 
















5900.406 




226 ( 6) 
0.810(5) 


















6082.899 




274 (21) 
0.812(12) 


















6108.403 






394 (35) 
0.844(14) 
















6133.957 








259 (65) 
0.789 (39) 














6210.663 


194 (24) 
0.252 (20) 


230 ( 9) 
0.813(6) 
















113 (27) 
0.767 (38) 


6248.274 




241 (7) 
0.819(5) 


397 (31) 
0.851 (12) 


343 (28) 
0.815(13) 












158 (13) 
0.782(13) 


6250.438 




306 ( 9) 
0.814(5) 


















6263.202 




255 (11) 
0.825 ( 7) 


















6287.20 




231 ( 8) 
0.823 ( 6) 




315(36) 
0.805 (19) 














6298.411 






399 (32) 
0.859 (13) 
















6341.492 






539 (54) 
0.871 (16) 














140 ( 9) 
0.790(11) 


6354.75 




241 (9) 
0.811 (6) 


















6365.540 




239 ( 6) 
0.814(4) 


619 (49) 
0.840(13) 


437 (38) 
0.865 (14) 














6465.218 


142 (18) 
0.345 (20) 


















127 (14) 
0.779 (18) 


6492.335 






489 (42) 
0.867 (13) 
















6495.576 




214(6) 
0.816(5) 


















6504.464 




240 ( 7) 
0.801 (5) 


















6514.959 


171 (23) 
0.335 (21) 




















6585.696 




191 (5) 
0.820 ( 5) 


















6588.020 




218(7) 
0.820 ( 6) 


















6637.187 


180(13) 
0.339(12) 




523 (43) 
0.852(13) 


432 (42) 
0.848 (16) 








527 (75) 
0.233 (24) 




174(16) 
0.798 (15) 


6637.96 


244(13) 


219(6) 


414 (28) 


484(19) 




1015 (104) 








111(6) 




0.346 ( 9) 


0.826 ( 5) 


0.821 (11) 


0.854 ( 6) 




0.826 ( 16) 








0.764 (10) 


6650.517 


241 (12) 


196 ( 4) 


480 (28) 


435 (21) 




865 ( 86) 








189 ( 7) 




0.358 ( 8) 


0.822 ( 3) 


0.848 ( 9) 


0.872( 8) 




0.811 ( 15) 








0.748 ( 7) 


6680.137 


157 (19) 
0.327(19) 




















6728.923 




214(6) 
0.819(5) 


















6740.078 




221 (7) 
0.841 (6) 


577 (52) 
0.821 (14) 


412 (26) 
0.872(10) 












221 (15) 
0.833 (11) 


6749.27 




262 (10) 
0.811 (6) 


















6764.596 




226 (10) 
0.822 ( 7) 


















6812.330 




208 ( 7) 
0.811 (5) 


















6816.004 




258 ( 8) 
0.813(5) 


484 (51) 
0.827 (17) 

















Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


Nd III 


3597.628 


137 ( 7) 






















0.484 ( 9) 




















o^ao aoo 
3603.983 


185 (15) 






















A JTl / I A \ 

KIA/2 (14) 




















3644.353 


148 ( 9) 






















A A O O /I At 

0.482 (10) 




















4570.660 














1013 (59) 






















0.6/6 ( 9) 








a n c a co 

4759.526 














119 (42) 






















0.712 ( 9) 








4769.622 














802 (51) 






















0.686 (10) 








4796.499 


212 ( 6) 






















A ^ ^ A / C\ 

0.439 ( 5) 




















4821.986 


209 ( 8) 












1039 (74) 










0.389 ( 6) 












0.641 (ll) 








491 1.653 














654 (33) 






















a n 1 o / o\ 

U. / lis ( s) 








4912.944 


241 ( 8) 






















0.429 ( 5) 




















4927.488 


208 ( 7) 






















0.477 ( 6) 




















5050.695 






375 (11) 




1058 (26) 




COO j"0 A \ 

583 (24) 






92(4) 








0.853 ( 5) 




1.157(4) 




0.751 ( 7) 






0.262 ( 8) 


5084.660 




203 ( 3) 


A C A / I i\\ 

454 (19) 




















A OO O / 0\ 

0.828 ( 2) 


0.852 ( 7) 
















5127.044 


163 ( 6) 




398 (10) 














86(2) 




0.498 (6) 




A O^A / A \ 

0.S60 ( 4) 














0.345 ( 5) 


5151.731 


215 (7) 


227 ( 3) 
















29(4) 




0.411 ( 5) 


A O C 1 / 0\ 

0.851 ( 2) 
















A O 1 A /A.A 

0.314 (26) 


coao ao i 

5203.924 


159 ( 6) 


251 (3) 


465 (13) 




A A O i A C \ 

948 (45) 




474 (35) 


ATA / A\ 

373 ( 9) 


A A 1 / I 0\ 

447 (18) 


128 ( 4) 




A CAO / f\ 

0.503 ( 6) 


O.o II ( 2) 


A OtZI 1 C\ 

0.86/ (5) 




1 ">0 1 / O s 

1.231 ( o) 




A O 1 A / 1 1\ 

O.ol9 (12) 


0.613 ( 4) 


A £00 / "7\ 

0.63o ( /) 


A OTO / C\ 

0.3 15 ( 5) 


5286.753 


267 ( 9) 


268 ( 5) 


A AO /^OA 

498 (23) 


541 (14) 


A 1 O 

918 (33) 


717 (37) 


O A C /O 0\ 

845 (38) 


AOO /'OC\ 

933 (25) 




28 ( 5) 




A 1 O 1 / H \ 

0.381 ( 6) 


A O A A 1 O \ 

0.844 ( 3) 


A O A O / 

0.843 ( 7) 


A A AO / A \ 

0.908 ( 4) 


1 AOA / £\ 

1 .039 ( 6) 


A O ~7 / 0\ 

0.867 ( 8) 


A /"I A / 

0.619 ( 7) 


A AA/ / A \ 

0.386 ( 4) 




0.718 (30) 


co a •! 1 i o 

5294.1 13 


188 ( 8) 


303 ( 4) 


(Ol / 1 A*\ 

432 (10) 


198 ( 8) 


OTO. /O A \ 

872 (34) 




C A 1 /A A\ 

541 (40) 


OAO / l\\ 

298 ( 9) 


^A1 /OA*\ 

601 (30) 


138 ( 5) 




A C Af\ f 1\ 

0.540 ( 7) 


0.888 ( 2) 


0.854 ( 4) 


1.100 ( 7) 


1.240 ( 6) 




0.864 (12) 


A O A / C\ 

0.630 ( 5) 


A 7MO / 0\ 

0.708 ( 8) 


0.340 ( 7) 


r I in AHA 

5410.099 


213 ( 7) 




O AO / 1 A\ 

392 (10) 




841 (15) 










47 ( 1) 




0.457 ( 5) 




0.847 ( 4) 




1.100 ( 3) 










A O O H t £\ 

0.286 ( 6) 


5566.015 




169 ( 3) 


363 (16) 


340 ( 9) 


781 (27) 










26 ( 2) 






A Q AC / O \ 

0.845 ( 3) 


0.864 ( 7) 


A A O IL i C\ 

0.926 ( 5) 


1 .052 ( 6) 










0.648 (17) 


5633.554 


241 ( 9) 


238 ( 4) 


iA^ /^A\ 

492 (22) 


AO£ i 1 A \ 

486 (14) 


1 A A A /O £I\ 

1049 (36) 


1 A "7 A / A £\ 

1070 (46) 








16(4) 




A A AC ( C\ 

0.405 ( 6) 


A OA<\ ( 1\ 

0.840 ( 3) 


A OCA 1 n\ 
0.059 ( /) 


A A O A 1 C\ 

0.930 (5) 


1.060 ( 6) 


A OAA / H\ 

0.Q90 ( /) 








A COO ( A \ 

0.5o2 (46) 


5677.179 


236 ( 9) 


164 ( 3) 


442 ( 1 6) 


421 (12) 


1190 (39) 


1160 (49) 


901 (48) 


565 (15) 


572 (40) 


52 ( 8) 




0.431 ( 6) 


A OCA / O \ 

0.854 ( 3) 


A O C O / £ \ 

0.852 ( 6) 


A AC A i C\ 

0.954 ( 5) 


1 AOO / C\ 

1.088 ( 5) 


A AO A / 

0.934 ( 7) 


A £.0£. / A\ 

0.686 ( 9) 


A A C A i A\ 

0.454 ( 4) 


0.590 (11) 


0.21 1 (26) 


5802.532 




220 ( 3) 


511 (19) 


441 (11) 


1137 (34) 


1055 (30) 


869 (45) 


503 (11) 


640 (44) 


66 ( 6) 






A O A "1 / O \ 

0.847 ( 3) 


0.847 ( 6) 


A ACO / A \ 

0.952 ( 4) 


1 A AO / C\ 

1.098 ( 5) 


A A 1 A S C\ 

0.914 ( 5) 


0.683 ( 8) 


0.454 ( 4) 


0.581 (11) 


0.297 (16) 


CO A C A") A 

5845.020 


22/ (8) 


275 ( 6) 


AOI t 1 0\ 

483 (18) 


369 (12) 


1 i at /oo\ 
1 10/ (32) 


1 ATO / A O \ 

10/2 (42) 








65 ( 5) 




0.466 ( 6) 


A O C C / A \ 

0.855 ( 4) 


0.850 ( 6) 


1 AAA / C\ 

1.009 ( 5) 


1.129 ( 6) 


A A O C / H\ 

0.935 ( 6) 








aoio f\ 

0.313 (12) 


5851.542 


261 ( 9) 


277 ( 8) 


461 (17) 


462 (14) 


1319 (49) 


1267 (55) 


863 (87) 


547 (18) 


663 (36) 


65 (8) 




0.431 ( 6) 


A O O ZT / C \ 

0.836 ( 5) 


0.851 ( 6) 


0.952 ( 5) 


1.107 ( 6) 


0.938 ( 7) 


0.679 (16) 


0.457 ( 5) 


A C C A / A\ 

0.554 ( 9) 


A O A O /O A\ 

0.243 (20) 


5987.683 


245 ( 8) 












/"AA i A A\ 

692 (44) 










A A C A 1 C\ 

0.464 ( 6) 












A TO A / 1 A\ 

0. /30 (10) 








6145.068 


194 ( 9) 


197 ( 7) 


266 (10) 


226 ( 8) 


1116 (38) 


OA A /TH 

829 (32) 




o o a 

324 (12) 


585 (43) 


109 ( 8) 




0.482 ( 8) 


0.857 ( 6) 


0.856 ( 6) 


1.015(6) 


1.144(5) 


0.982 ( 6) 




0.518 (6) 


0.624(12) 


0.270 ( 1) 


6327 265 








311 (10) 




999 (56) 




375 (15) 


543 (S51 


94 ( 7) 










1.030(5) 




0.971 (9) 




0.514(6) 


0.628 (16) 


0.248 (13) 


6550.231 








371 (12) 


1310(50) 








657 (60) 


102 ( 6) 










1.036(5) 


1.154(6) 








0.646(15) 


0.298 ( 9) 


6690.830 


235 (11) 


227 ( 7) 


481 (32) 


469 (21) 


1209 (96) 


1295 (70) 




761 (33) 


445 (139) 


40(8) 




0.399 ( 8) 


0.834 ( 5) 


0.848 (11) 


0.921 ( 7) 


1.056(13) 


0.879 ( 9) 




0.391 (7) 


0.550 (49) 


0.751 (34) 


Eu II 


5818.746 










515 (50) 


482 (45) 




















0.940 (16) 


0.721 (15) 










5872.978 




136(6) 






















0.786 ( 8) 


















6049.513 


74(12) 






















0.243 (26) 




















6173.029 


79(10) 


75(5) 




82(12) 




364 (30) 












0.248 (21) 


0.800(11) 




0.570 (25) 




0.736(13) 










6303.423 






84 (26) 




555 (46) 


356 (52) 




145 (21) 


282(71) 










0.788 (50) 




0.950 (13) 


0.740 (24) 




0.164(23) 


0.423 (40) 




6437.640 


109(10) 


57 (41) 


109(17) 


117(12) 


662 (41) 


469 (34) 






525 (77) 


151(5) 




0.247(14) 


0.807(11) 


0.800 (26) 


0.615 (17) 


0.943 (10) 


0.737 (12) 






0.397 (23) 


0.009 ( 6) 


6645.110 


110(8) 


38(4) 


105(13) 


172(9) 


600 (64) 


489 (48) 




227 ( 9) 


300 (58) 


191 (6) 




0.252(12) 


0.800(18) 


0.858 (20) 


0.657 ( 8) 


0.986(17) 


0.723 (16) 




0.184(6) 


0.472 (31) 


-0.001 (5) 


Tbll 


6194.488 




250 (56) 






















0.956 (35) 



















Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 


Tb III 


5505.409 


326 (19) 




1129 (30) 


501 (22) 


1242 (31) 






985 (17) 


534 (28) 


14(2) 




0.884 ( 9) 




0.957 ( 4) 


1.404 ( 7) 


1.228 ( 4) 






0.864 ( 3) 


0.675 ( 9) 


0.543 (31) 


5847.232 


359 (23) 




988 (33) 


519 (30) 








1150(45) 


474 (89) 


103 ( 9) 




0.869 (10) 




0.926 ( 5) 


1.387 ( 9) 








0.857 ( 6) 


0.652 (30) 


0.007 (15) 


6092.897 


218(13) 
0.841 (10) 




991 (34) 
0.963 ( 6) 




677 (35) 
1.215 ( 8) 










140 (10) 
0.078 (12) 


6323.620 




654 (17) 


1141 (51) 


556 (45) 


1 124 (77) 


507 (56) 




965 (51) 


495(107) 


202 (19) 






1.119 ( 4) 


0.942 ( 7) 


1.401 (13) 


1.237 (11) 


1.109 (18) 




0.826 ( 9) 


0.799 (34) 


t\ AO O / \ C\ 

0.038 (15) 


6511.060 


329 (23) 
0.824 (1 1) 


480 (11) 
1.058 ( 4) 






833 (75) 
1.247 (14) 


459 (45) 
1.134 (16) 






472(117) 
0.712 (41) 




6511.788 




551 (13) 
1.033 ( 4) 






854 (67) 
1.194 (13) 


449 (57) 
1.112 (20) 










6591.215 






1022 (65) 
0.977 (10) 
















6687 698 

/ .U7U 


197 (23) 
0.846(18) 




913 (551 
0.936 ( 9) 




1 067 (?271 

1 \}\.J I 1 J 

1.209(34) 


616 (1461 
1.105 (38) 








126 (18) 
-0.031 (23) 


6716.566 






913 (42) 
0.932 ( 7) 










1041 (175) 
0.818 (26) 






6832.056 




544 (19) 
1.042(6) 


















Dy II 


3407.796 


230(14) 
0.485 (10) 




















3429.443 


382 (24) 
0.406(10) 




















3434.369 


283 (14) 
0.449 ( 8) 




















3534.960 


227 (13) 
0.470 ( 9) 




















3538.519 


292(14) 
0.467 ( 8) 




















3550.218 


315(19) 
0.453 (10) 




















3563.148 


317 (18) 
0.443 ( 9) 




















3602.823 


178 (39) 
0.315(35) 




















4923.167 


201 (13) 
0.372(11) 




















5090.386 






327 (23) 
0.826(11) 
















5164.115 








433 (31) 
0.825 (11) 












144(12) 
0.704(14) 


5169.688 






356 (20) 
0.843 ( 9) 


368 (14) 
0.841 ( 6) 












84(4) 
0.785 ( 9) 


5272.249 




129 ( 3) 
0.811 (4) 


















5279.697 




174 ( 4) 
0.816(4) 


















5284.989 




127 ( 3) 
0.761 (5) 


















5337.427 




145 ( 6) 
0.818(7) 


















5399.936 




125 ( 6) 
0.829 ( 8) 


















5426.695 






335 (45) 
0.811 (21) 


473 (25) 
0.816(9) 




268 (66) 
0.536 (39) 










5443.343 




159 ( 6) 
0.802 ( 7) 


















5471.909 




150 ( 4) 
0.815 (5) 








214(99) 
0.629 (70) 










5515.405 




216(5) 
0.807 ( 4) 


















5641.504 




121 (6) 
0.816(8) 


















5868.107 






419 (57) 
0.831 (22) 














77 (14) 
0.867 (30) 


5924.558 




120(7) 
0.837 (10) 


















6196.234 


160(31) 
0.329(31) 


230 ( 8) 
0.798 ( 6) 


















6654.228 




161 (9) 
0.830 ( 9) 
















124(16) 
0.837 (22) 


6700.639 


193 (30) 
0.320 (24) 




426 (68) 
0.869 (25) 


311 (52) 
0.836 (27) 












171 (23) 
0.803 (22) 


Dy III 


3919.398 


278 (14) 
0.490 ( 8) 




















5730.284 


255 (12) 


129 ( 3) 


403 (23) 


496(15) 


683 (42) 


705 (61) 


356 (61) 


856 (42) 


461 (76) 


39(4) 



Table 3. Continued. 



Wavelength 


Amplitude (m s 1 ) / Phase (period fraction) 


(A) 


HD 24712 


HD 101065 


HD 122970 


HD 134214 


HD 12932 


HD 19918 


HD 201601 


HD 128898 


HD 9289 


HD 137949 




0.451 (8) 


0.838 (4) 


0.848 ( 9) 


0.893 (5) 


1.068(10) 


0.866 (14) 


0.671 (28) 


0.413 ( 8) 


0.518(27) 


0.446(18) 


Th III 


5376.130 
6599.482 


39 (22) 
42(19) 


388 (15) 
1.120(6) 

224 (13) 
1.086(9) 


670 (64) 
1.088 (15) 
807 (120) 
1.047 (23) 


273 (17) 
1.425 (11) 

402 (53) 
1.437 (21) 












186(13) 
-0.111 (7) 

302(12) 
0.001 ( 7) 



